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In a recent paper [l[ it was shown that in fluxless M theory vacua with at least two hidden 
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sectors undergoing strong gauge dynamics and a particular form of the Kahler potential, all 
moduli are stabilized by the effective potential and a stable hierarchy is generated, consis- 
tent with standard gauge unification. This paper explains the results of [l| in more detail 
and generalizes them, finding an essentially unique de Sitter (dS) vacuum under reasonable 
conditions. One of the main phenomenological consequences is a prediction which emerges 
from this entire class of vacua: namely gaugino masses are significantly suppressed relative 
to the gravitino mass. We also present evidence that, for those vacua in which the vacuum 
energy is small, the gravitino mass, which sets all the superpartner masses, is automatically 
in the TeV - 100 TeV range. 
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I. INTRODUCTION AND SUMMARY 

There are many good reasons why we study string theory as a theory of particle physics. One 
of these, discovered some twenty or so years ago [2], is that a simple question, "what properties do 
four-dimensional heterotic string vacua generically have?" has an extremely compelling answer: non- 
Abelian gauge symmetry, chiral fermions, hierarchical Yukawa couplings and dynamical supersymmetry 
breaking. One can also add other important properties such as gauge coupling unification and doublet- 
triplet splitting. Furthermore, after the dust of the string duality revolution settled, a similar picture was 
discovered in other perturbative corners of the landscape, eg Type IIA, Type IIB and M theory. The 
above four properties are the most important properties of the Standard Model. The fifth one - gauge 
coupling unification - is an important feature of the MSSM and low energy supersymmetry. 

There is another crucial feature of the Standard Model: namely its overall mass scale, which is of 
order M\y as opposed to some other scale such as Mp or Mqut- This property of the Standard Model 
is much less well understood in string theory or otherwise, and will be the focus of this paper. 

In string/M theory, masses and coupling constants, including Myy are all functions of the moduli 
field vevs. Thus the hierarchy problem in string theory is double edged: one has to both stabilize all 
the moduli and generate the hierarchy simultaneously. In recent years there has been progress in moduli 
stabilization via fluxes and other effects (for a recent review see p) ) , and with the imminent arrival of the 
LHC it is appropriate to address the hierarchy problem in this context. After all, any given string/M 
theory vacuum either will or will not be consistent with the LHC signal, but one cannot even begin to 
address this in a meaningful way if the hierarchy is not understood. 

In Type IIB vacua, moduli are stabilized via a combination of fluxes and quantum corrections HQ. In 
these vacua, hierarchies arise in three ways: warp factors 0,3 as in 0, the presence of non-perturbative 
effects or by fine tuning the large number of fluxes. However, flux vacua in Type IIA M theory 
0] and heterotic string theory have the property that the (currently known and understood) fluxes 
roughly are equal in number to the number of moduli. This leads to a large value of the superpotential, 
and consequently, if the volume of the extra dimensions is not huge, to a large gravitino mass. This 
tends to give a large mass to all scalars via the effective 4d supergravity potential, which leads to either 
M\y = or some other large value such as Mp or Mqut- Therefore, for these vacua we require a good 
idea for generating and stabilizing the hierarchy. 

Thus far, there has been essentially one good idea proposed to explain the relatively small value of the 
weak scale. This is that the weak scale might be identified with, or related to, the strong coupling scale 
of an asymptotically free theory which becomes strongly coupled at low energies and exhibits a mass gap 
at that strong coupling scale. Holographically dual to this is the idea of warped extra dimensions [7]. 
Strong dynamics (or its dual) can certainly generate a small scale in a natural manner, but can it also 
be compatible with the stabilization of all the moduli fields? 

One context for this question, which we will see is particularly natural, is M theory compactification 
on manifolds X of GVholonomy without fluxes. In these vacua, the only moduli one has are zero modes of 
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the metric on X, whose bosonic superpartners are axions. Thus each moduli supermultiplet has a Peccei- 
Quinn shift symmetry (which originates from 3-form gauge transformations in the bulk lid super gravity). 
Since such symmetries can only be broken by non-perturbative effects, the entire moduli superpotential 
W is non-perturbative. In general W can depend on all the moduli. Therefore, in addition to the small 
scale generated by the strong dynamics we might expect that all the moduli are actually stabilized. This 
paper will demonstrate in detail that this is indeed the case. 

Having established that the basic idea works well, the next question we address is "what are the 
phenomenological implications?" Since string/M theory has many vacua, it would be extremely useful 
if we could obtain a general prediction from all vacua or at least some well-defined subset of vacua. 
Remarkably, we are able to give such a prediction for all fluxless M theory vacua within the supergravity 
approximation 1 with at least two hidden sectors undergoing strong gauge dynamics and a particular form 
of the Kahler potential as in (pQ): gaugino masses are generically suppressed relative to the gravitino mass. 

A slightly more detailed elucidation of this result is that in all de Sitter vacua within this class, gaugino 
masses are always suppressed. In AdS vacua - which are obviously less interesting phenomenologically - 
the gaugino masses are suppressed in 'most' of the vacua. This will be explained in more detail later. 

The reason why we are able to draw such a generic conclusion is the following: any given non- 
perturbative contribution to the superpotential depends on various constants which are determined by 
a specific choice of GVmanifold X. These constants determine entirely the moduli potential. They are 
given by the constants bk which are related to the one-loop beta-function coefficients, the normalization 
of each term and the constants a, (see ([2])) which characterize the Kahler potential for the moduli. 
Finally, there is a dependence on the gauge kinetic function, and in M theory this is determined by a set 
of integers N{ which specify the homology class of the 3-cycle on which the non-Abelian gauge group is 
localized. Rather than study a particular X which fixes a particular choice for these constants we have 
studied the effective potential as a function of the (Ak, 6fc, Qj, iVj) . The result of gaugino mass suppression 
holds essentially for arbitrary values of the (Ak,bk,a,i,Ni), at least in the supergravity regime where we 
have been able to calculate. Thus, any Cr2-manifold which has hidden sectors with strong gauge dynamics 
will lead to suppressed gaugino masses. 

At a deeper level, however, the reason that this works is that the idea of strong gauge dynamics to 
solve the hierarchy problem is a good and simple idea which guides us to the answers directly. If one's 
theory does not provide a simple mechanism for how the hierarchy is generated, then it is difficult to see 
how one could obtain a reliable prediction for, say, the spectrum of beyond the Standard Model particles. 
In a particular subset of Type IIB compactifications, Conlon and Quevedo have also discovered some 
general results [llj. In fact, they remarkably also find that gaugino masses are suppressed at tree level, 
though the nature of the suppression is not quite the same. Some heterotic compactifications also exhibit 
a suppression of tree- level gaugino masses [12|] . 

The suppression of gaugino masses relative to applies for all vacua in the supergravity regime 

arising out of these compactifications, independent of the value of However, in a generic vacuum 

the cosmological constant is too large. If we therefore consider only those vacua in which the cosmological 
constant is acceptable at leading order, this constrains the scale of m 3 /2 further. Remarkably, we find 
evidence that for such vacua, 7/13/2 is of order 1 — 100 TeV. This result certainly deserves much further 
investigation. 

The fact that such general results emerge from these studies makes the task of predicting implications 
for various collider observables as well as distinguishing among different vacua with data from the LHC 
(or any other experiment) easier. A more detailed study of the collider physics and other phenomenology 
will appear in the future [13j |. However, as we will see in section IVIIIl it could be quite easy to distinguish 



1 to which we are restricted for calculability 
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Type IIB and M theory vacua using the forthcoming LHC data. 

This paper is the somewhat longer companion paper to [l|. Given its length we thought that it would 
be worthwhile to end this introduction with a guide to its contents. Much of the bulk of the paper 
is devoted to analyzing and explaining the details of why the potential generated by strong dynamics 
in the hidden sector has vacua in which all moduli are stabilized. At first we begin with the simplest 
non-trivial example, two hidden sectors, without any charged matter: only gauge bosons and gauginos. 
Section [TT] calculates the moduli potential in this case. Section IIIII analyzes its supersymmetric vacua: 
these are all isolated with a negative vacuum energy. Section IIVI describes explicit examples realizing 
the vacuum structure of Sections IIIII and [V] Section [V] describes the vacua which spontaneously break 
supersymmetry. These also have negative vacuum energy and all moduli stabilized. Section IVII goes 
on to consider more complicated hidden sectors, where it is argued that metastable de Sitter vacua can 
also occur under very reasonable conditions, and that the metastable de Sitter vacuum obtained for a 
given GVmanifold is essentially unique. In section IVlII we study the distribution of 7713/2- In particular 
for the de Sitter vacua it is shown that requiring the absence of a large cosmological constant fixes the 
gravitino mass to be of 0(1-100) TeV. Section [Villi discusses phenomenology. In particular, we explain 
the suppression of gaugino masses and discuss the other soft SUSY breaking couplings. We conclude in 
section IIXI followed by an Appendix which discusses the Kahler metric for visible charged matter fields 
in M theory. 



II. THE MODULI POTENTIAL 

In this section we quickly summarize the basic relevant features of G2-compactifications, setup the 
notation and calculate the potential for the moduli generated by strong hidden sector gauge dynamics. 

In M theory compactifications on a manifold X of GVholonomy the moduli are in correspondence 
with the harmonic 3-forms. Since there are N = 63(A) such independent 3-forms there are N moduli 
Zi = ti + isi . The real parts of these moduli ti are axion fields which originate from the 3-form field C in 
M theory and the imaginary parts s, are zero modes of the metric on X and characterize the size and 
shape of X. Roughly speaking, one can think of the Sj's as measuring the volumes of a basis of the N 
independent three dimensional cycles in X. 

Non-Abelian gauge fields are localized on three dimensional submanifolds Q of X along which there is 



an orbifold singularity 1J] while chiral fermions are localized at point-like conical singularities 15|, [l6|, [17|] . 
Thus these provide M theory realizations of theories with localized matter. A particle localized at a point 
p will be charged under a gauge field supported on Q if p E Q. Since generically, two three dimensional 
submanifolds do not intersect in a seven dimensional space, there will be no light matter fields charged 
under both the standard model gauge group and any hidden sector gauge group. Supersymmetry breaking 
is therefore gravity mediated in these vacua. 2 

In general the Kahler potentials for the moduli are difficult to determine in these vacua. However a 
set of Kahler potentials, consistent with GVholonomy and known to describe accurately some explicit 



examples of G2 moduli dynamics were given in 18[]. These models are given by 



K = -31n(4vr 1/3 y x ) (1) 



2 This is an example of the sort of general result one is aiming for in string/Af theory. We can contrast this result with 
Type IIA vacua. Here the non-Abelian gauge fields are again localized on 3-cycles, but since generically a pair of three 
cycles intersect at points in six extra dimensions, In Type IIA supersymmetry breaking will generically be gauge mediated. 
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where the volume in 11-dimensional units as a function of s, is 

N N 

V x = Y[s?, with ^ = 7/3 (2) 

i=i i=i 

We will assume that this iV-parameter family of Kahler potentials represents well the moduli dynamics. 
More general Kahler potentials outside this class have the volume functional multiplied by a function 
invariant under rescaling of the metric. It would be extremely interesting to investigate the extension of 
our results to these cases. 

As motivated in the introduction, we are interested in studying moduli stabilization induced via strong 
gauge dynamics. We will begin by considering hidden sector gauge groups with no chiral matter. Later 
sections will describe the cases with hidden sector chiral matter. 

In this 'no matter' case, a superpotential (in units of m? p ) of the following form is generated 

M 

W = Y^ A k e ibkf * (3) 
k=i 

where M is the number of hidden sectors undergoing gaugino condensation, bk = t 2 - with Ck being the 
dual coxeter numbers of the hidden sector gauge groups, and are numerical constants. The A^ are 
RG-scheme dependent and also depend upon the threshold corrections to the gauge couplings; the work 
of [3] shows that their ratios (which should be scheme independent) can in fact take a reasonably wide 
range of values in the space of M theory vacua. We will only consider the ratios to vary from O(0. 1-10) 
in what follows. 

The gauge coupling functions fk for these singularities are integer linear combinations of the z% , because 
a 3-cycle Q along which a given non-Abelian gauge field is localized is a supersymmetric cycle, whose 
volume is linear in the moduli. 

N 

/* = J>jV (4) 
i=i 

Notice that, given a particular GVmanifold X for the extra dimensions, the constants (a^, bk, A^, Nf) are 
determined. Then, the Kahler potential and superpotential for that particular X are completely deter- 
mined by the constants (ctj, bk,Ak, N*). This is as it should be, since M theory has no free dimensionless 
parameters. 

We are ultimately aiming for an answer to the question, "do M theory vacua in general make a 
prediction for the beyond the standard model spectrum?". For this reason, since a fluxless M theory 
vacuum is completely specified by the constants (ai,bk, A^, Nf) we will try as much as possible not to 
pick a particular value for the constants and try to first evaluate whether or not there is a prediction for 
general values of the constants. Our results will show that at least within the supergravity approximation 
there is indeed a general prediction: the suppression of gaugino masses relative to the gravitino mass. 

At this point the simplest possibility would be to consider a single hidden sector gauge group. Whilst 
this does in fact stabilize all the moduli, it is a) non-generic and b) fixes the moduli in a place which is 
strictly beyond the supergravity approximation. Therefore we will begin, for simplicity, by considering 
two such hidden sectors, which is more representative of a typical G2 compactification as well as being 
tractable enough to analyze. The superpotential therefore has the following form 



W np = A ie iblh +A 2 e ib2h 



(5) 



7 



The metric corresponding to the Kahler potential (pQ) is given by 



K i3 = J^ S i3- ( 6 ) 



The ]\f = 1 supergravity scalar potential given by 

V = e K {K rj FiFj -3\W\ 2 ), (7) 

where 

Fi = diW + (diK)W , (8) 
can now be computed. The full expression for the scalar potential is given by 

2 N 2 

v = i^EE^K* +3 ) 6 ^^^^ 

X k=l i=l k=l 
N 2 

+2cos[(6 1 iV 1 - b 2 N 2 ) ■ t}J2^i II "ihA k e- bki;k - s (9) 

8=1 k=l 

+3cos[(6iiV 1 - b 2 N 2 ) • t] U + Y, b k " k • a J JJ V^'l 

V fc=i / j=i 



where we introduced a variable 



such that 

Im/ fc = z7 fc -a. (11) 

By extremizing Q with respect to the axions ti we obtain an equation 

sin^N 1 -b 2 N 2 ) -t] = 0, (12) 

which fixes only one linear combination of the axions. In this case 

cos[(6iiV 1 - b 2 N 2 ) -t] = ±1. (13) 

It turns out that in order for the potential @ to have minima, the axions must take on the values such 
that cos[(6iiV 1 — b 2 N 2 ) ■ t ] = —1 for Ai,A 2 > 0. Otherwise the potential has a runaway behavior. After 
choosing the minus sign, the potential takes the form 

1 2 N 2 

V = I^W EE*** 0* + 3) hAle-^ Ks + 3J2^ 2bki;Ka (14) 



* k=l t=l k=l 



N 2 / 2 \ 2 

2 E °< II ^Mfce" 6 ***' " - 3 2 + ^ 6 fc z? fc • a [J 
i=i fc=i V fe=i / i=i 



In the next section we will go on to analyze the vacua of this potential with unbroken super symmetry. 
The vacua in which supersymmetry is spontaneously broken are described in sections V and VI. 
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III. SUPERSYMMETRIC VACUA 



In this section we will discuss the existence and properties of the supersymmetric vacua in our theory. 
This is comparatively easy to do since such vacua can be obtained by imposing the supersymmetry 
conditions instead of extremizing the full scalar potential (|14p . Therefore, we will study this case with 
the most detail. Experience has also taught us that potentials possessing rigid isolated supersymmetric 
vacua, also typically have other non-super symmetric vacua with many qualitatively similar features. 

The conditions for a supersymmetric vacuum are: 

Fi = 0, 

which implies 

{b x N} + p-)A, + (b 2 N 2 + p-)A 2 [cosfM 1 - b 2 N 2 ) • t] e^ 1 -^ 2 )- 

+ ism[{b 1 N 1 -b 2 N 2 )-t]e ibl ^- b2 ^ 2) - g ] = (15) 

Equating the imaginary part of (|15|) to zero, one finds that 

sinK&iiV 1 - b 2 N 2 ) ■ t] = 0, 
which implies cos[(6iiV 1 - b 2 N 2 ) ■ t] = ±1 . (16) 

For A\,A 2 > 0, a solution with positive values for the moduli (sj) exists when the axions take on the 
values such that cos[(6iiV 1 — b 2 N 2 ) ■ t] = —1 . Now, equating the real part of (fT5|) to zero, one obtains 

(hNl + p-)A 1 + (b 2 N 2 + p-)A 2 ^-^K = . (17) 

This is a system of N transcendental equations with N unknowns. As such, it can only be solved 
numerically, in which case, it is harder to get a good understanding of the nature of solutions obtained. 
Rather than doing a brute force numerical analysis of the system (|17p it is very convenient to introduce 
a new auxiliary variable a to recast (fTTl) into a system of linear equations with iV unknowns coupled to 
a single transcendental constraint as follows: 

a(b 1 Nl + ^)-(b 2 Nl + ^) = (18) 
A\ a 

The system of linear equations (I18p can then be formally solved for in terms of {b\,b 2 , Nl,Nf, a,} and 



a: 



3ai(a - 1} • i = l,2,..,N. (20) 



2 {hWa - b 2 m 



One can then substitute the solutions for Sj into the constraint (|19p and self-consistently solve for the 
parameter a in terms of the input quantities {A\, A 2 ,b±,b 2 , N^, N 2 ,a>i}. This, of course, has to be done 
numerically, but we have indeed verified that solutions exist. Thus, we have shown explicitly that the 
moduli can be stabilized. We now go on to discuss the solutions, in particular those which lie within the 
supergravity approximation. 
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A. Solutions and the Supergravity Approximation 



Not all choices of the constants {A±, A 2 , 61, b 2 , N^, Nf, a{\ lead to solutions consistent with the ap- 
proximation that in the bulk of spacetime, eleven dimensional supergravity is valid. Although this is not 
a precisely (in the numerical sense) defined approximation, a reasonable requirement would seem to be 
that the values of the stabilized moduli (sj) obtained from (|20p are greater than 1. It is an interesting 
question, certainly worthy of further study, whether or not this is the correct criterion. In any case, this 
is the criterion that we will use and discuss further. 

From ()20p and (|19p . and requiring the Sj to be greater than 1, we get the following two branches of 
conditions on parameter a : 

a) ->1; rmn{^ = l,N}>a> max + ^ = l,N } 

b) £<1; ma X {^§ ;i= T^} < a < mm ffffi + f ^ ;i = T^} (21) 
Ax biNl bxNl + 3aj/2 

Notice that the solution for Sj (|20p has a singularity at a = iNct- This can be seen clearly from Figure 

1 i 

dl]). We see that the modulus s% (> 0) falls very rapidly as one moves away from the vertical asymptote 
representing the singularity and can become smaller than one very quickly, where the supergravity ap- 
proximation fails to be valid. The relative location of the singularities for different moduli will turn out 




FIG. 1: Positive values of Si plotted as a function of a for a case with two condensates and three bulk moduli for the 
following choice of constants 61 = |f ,62 = fg-, iV/ = {1,2,2}, AT? = {2,3, 5}, a; = {1,1/7,25/21}. The qualitative feature 
of this plot remains the same for different choices of constants as well as for different i. The vertical line is the locus for 

a = b . 2l ^ T \ , where the denominator of (120 11 vanishes. 

01 Iyj~ ' 

to be very important as we will see shortly. From (|21h . we know that there are two branches for allowed 
values of a. Here we consider branch a) for concreteness, branch b) can be analyzed similarly. Figure (|2|) 
shows plots for A-ijA\ and s« as functions of a for a case with two condensates and three bulk moduli. 
The plots are for a given choice of the constants {b\,b2,Nl ,Nf,ai, i = 1,2,3}. The qualitative feature 
of the plots remains the same even if one has a different value for the constants. 
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FIG. 2: Left - A2/A1 plotted as a function of a for a case with two condensates and three bulk moduli. The function 
diverges as it approaches the loci of singularities of 1|20[) . viz. a = . 

Right - Positive Si, i = 1, 2, 3 for the same case plotted as functions of a. si is represented by the solid curve, S2 by the long 
dashed curve and S3 by the short dashed curve. The vertical lines again represent the loci of singularities of (|20[l which the 
respective moduli Si asymptote to. The horizontal solid (red) line shows the value unity for the moduli, below which the 
supergravity approximation is not valid. 

Both plots are for 61 = §|, 62 = ff , Nl = {1, 2, 2}, iVf = {2, 3, 5}, a* = {1, 1/7, 25/21}. 



Since the Si fall very rapidly as one goes to the left of the vertical asymptotes, there is a small 
region of a between the origin and the leftmost vertical asymptote which yields allowed values for all 
Si > 1. Thus, for a solution in the supergravity regime all (three) vertical lines representing the loci of 
singularities of the (three) moduli Si should be (sufficiently) close to each other. This means that the 

positions of the vertical line for the ith modulus (a = ) and the jth modulus (a = ) can not be 

too far apart. This in turn implies that the ratio of integer coefficients (Nf /Nf) and (Nj /Nj) for the ith 
and jth modulus cannot be too different from each other in order to remain within the approximation. 
Effectively, this means that the integer combinations in the gauge kinetic functions ([!]) of the two hidden 
sector gauge groups in © can not be too linearly independent. We will give explicit examples of G2 
manifolds in which {N} /Nf) and (Nj/Nf) are the same for all i and j, so the constraint of being within 
the supergravity approximation is satisfied. 

We now turn to the effect of the other constants on the nature of solutions obtained. From the top 
right plot in Figure ([3j), we see that increasing the ranks of the gauge groups while keeping them close 
to each other (with all other constants fixed) increases the size of the moduli in general. On the other 
hand, from the bottom left plot we see that introducing a large difference in the ranks leads to a decrease 
in the size of the moduli in general. Hence, typically it is easier to find solutions with comparatively 
large rank gauge groups which are close to each other. The bottom right plot shows the sizes of the 
moduli as functions of A2/A1 while keeping the ranks of the gauge groups same as in the top left plot 
but changing the integer coefficients. We typically find that if the integer coefficients are such that the 
two gauge kinetic functions are almost dependent, then it is easier to find solutions with values of moduli 
in the supergravity regime. 

The above analysis performed for three moduli can be easily extended to include many more moduli. 
Typically, as the number of moduli grows, the values of at in (|20p decrease because of ([2]). Therefore the 
ranks of the gauge groups should be increased in order to remain in the supergravity regime as one can 
see from the structure of (|20p . At the same time, for reasons described above, the integer combinations 
for the two gauge kinetic functions should not be too linearly independent. In addition, the integers Nf 
should not be too large as they also decrease the moduli sizes in (|20p . 
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Si,S 2 ,S 3 Si,S 2 ,S 3 




FIG. 3: Plots of positive Si, i = 1, 2, 3 as functions of A2/A1. 

Top Left: Same choice of constants as in Figure©, i.e. 61 = §|, hi = §|, iVj = {1, 2, 2}, AT? = {2, 3, 5}, <Zj = {1, 1/7, 25/21}. 
Top Right: We increase the ranks of the gauge groups but keep them close (keeping everything else same) - 61 = fe ) fr 2 = 2s. 
Bottom Left: We introduce a large difference in the ranks of the gauge groups (with everything else same) - 61 = fe, 62 = §§■■ 
Bottom Right: We keep the ranks of the gauge groups as in Top Left but change the integer coefficients to = 
{1,2,2}, N} = {3,3,4}. 



What happens if some of the integers A^ 1 or Nf are zero. Figure H] corresponds to this type of a 
situation when the integer combinations are given by A^ 1 = {1, 0, 1}, N± ={1,1,1}. As we can see from 
the plots, all the moduli can still be stabilized although one of the moduli, namely S2 is stabilized at 
values less than one in 11-dim Planck units. This gets us back to the previous discussion as to when the 
supergravity approximation can be valid. We will not have too much to say about this point, except to 
note that a) the volume of X can still be large (([2]) is large, greater than one in 11-dim Plank units), 
b) the volumes of the associative three-cycles Qk which appear in the gauge kinetic function (JH), i.e. 
Vol(Qk) = Yli^L^i s i can a ^ so be large and c) that the top Yukawa in these models comes from a small 
modulus vev [13] . From Figure 0] we see that although the modulus S2 is always much smaller than 
one, the overall volume of the manifold Vx represented by the solid red curve is much greater than one. 
Likewise, the volumes of the associative three cycles Vol(Qi) = s\ + S3 and Vol{Q2) = s\ + S2 + S3 are 
also large. Therefore if one interprets the SUGRA approximation in this way, it seems possible to have 
zero entries in the gauge kinetic functions for some of the moduli and still stabilize all the moduli, as 
demonstrated by the explicit example given above. In general, however, there is no reason why any of 
the integers should vanish in the basis in which the Kahler metric is given by ©. 



12 




FIG. 4: Plots of positive Si,i = 1,2,3 as functions of Aa/Ai. The constants are 6j = |jj,fe 2 = ff.iV? = {1,0, l},iV? = 
{1, 1, 1}, a» = {1/10, 1, 37/30}. Si is represented by the solid curve, S2 by the long dashed curve and S3 by the short dashed 
curve. The red curve represents the volume of the internal manifold as a function of A2 /A\ . 
Right - the same plot with the vertical plot range decreased. 



B. Special Case 



A very interesting special case arises when the gauge kinetic functions /1 and fi in ([5]) are equal. 
Since in this case N\ = Nf, the moduli vevs are larger in the supersymmetric vacuum; hence this case 
is representative of the vacua to be found within the supergravity approximation. Even though this is a 
special case, in section HVl we will describe explicit examples of G2 manifolds in which = Nf. 

In the special case, we have 

Nl = Nf = N, , (22) 

and therefore 

12 _ NiSi 

Vi =Vi =v i = (23) 

For this special case, the system of equations ()1T|) can be simplified even further. We have 

{bxui + ^)At - {b 2 Vi + ^Ate^-W = (24) 

with Vi actually independent of i. Thus, we are left with just one simple algebraic equation and one 
transcendental constraint. The solution for is given by : 

Ui = t : ~^r~h — rr > ( 25 ) 

2(ab 1 - 6 2 ) 



with 



A\ a 

Since z/, is independent of i, it is also independent of the number of moduli N. In Figure ([5]) we plotted 
v as a function of A2/A1 when the hidden sector gauge groups are SU(5) and 5Z7(4). Notice that here 
the ranks of the gauge groups don't have to be large for the moduli to be greater than one. This is in 
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contrast with the linearly independent cases plotted in Figure ([3|). Once v is determined in terms of 
A2/A1, the moduli are given by: 

Si = ^. (27) 

Therefore, the hierarchy between the moduli sizes is completely determined by the ratios Oj/iVi for 
different values of %. In addition, from Figure [5] it can be seen that v keeps increasing indefinitely if we 




FIG. 5: Plot of v as a function of A2/A1 for the choice 61 = 62 = The solid (red) curve represents the exact 
numerical solution whereas the black dashed curve is the leading order approximation given by 



keep increasing AijA\ (though theoretically there may be a reasonable upper limit for A2/A1), which is 
not possible for the general case as there are N u^s. This implies that it is possible to have a wide range 
of the constants which yield a solution in the supergravity regime. 

Although the numerical solutions to the system (|25ti26|) described above are easy to generate, having 
an explicit analytic solution, even an approximate one, which could capture the dependence of v on the 
constants A2/A1, b\ and 62 would be very useful. 

Fortunately there exists a good approximation, namely a large v limit, which allows us to find an ana- 
lytical solution for v in a straightforward way. Expressing a from (|25p , in the leading order approximation 
when v is large we obtain 

a(°) = ^. (28) 
01 

After substituting (|28p into (|26p we obtain the approximate solution for v in the leading order: 

where the last expression corresponds to SU(P) and SU{Q) hidden sector gauge groups. For the moduli 
to be positive either of the two following conditions have to be satisfied 



a) AiQ < A 2 P; P>Q 

b) A 1 Q>A 2 P; P<Q. 



(30) 
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From the plots in Figure[5]we notice that the above approximation is fairly accurate even when v is 0(1). 
This is very helpful and can be seen once we compute the first subleading contribution. By substituting 
(|29[) back into (|25|) and solving for a we now have up to the first subleading order: 

a = a<°> + a™ = f F ~ V (31) 

It is then straightforward to compute v which includes the first subleading order contribution 

y = y m + ,w ± p r Q ,. (32) 

\AiQ 

We can now examine the accuracy of the leading order approximation when v is O(l) by considering 
the region where the ratio A<ij A\ is small. A quick check for the SU(5) and SU(4) hidden sector gauge 
groups chosen in the case presented in Figure [5] yields for A^jA\ = 4: 

a = a (o) +0,(4) = 1.25-0.136, (33) 
i/ = i/(°) = 2.195 - 0.148. (34) 

which results in a 12% and 7% error for cr ) and respectively. The errors get highly suppressed when 
v becomes 0(10) and larger. Also, when the ranks of the gauge groups SU(P) and SU(Q) are 0(10) 
and P — Q is small, the ratio A2/A1 can be O(l) and still yield a large v. The dependence of v on the 
constants in (|29p is very similar to the moduli dependence obtained for SUSY Minkowski vacua in the 
Type IIB racetrack models |2o| . 

We have demonstrated that there exist isolated supersymmetric vacua in M theory compactifications 
on G2-manifolds with two strongly coupled hidden sectors which give non-perturbative contributions to 
the superpotential. Given the existence of supersymmetric vacua, it is very likely that the potential 



also contains non-supersymmetric critical points. Previous examples have certainly illustrated this [181 ]. 
Before analyzing the non-supersymmetric critical points, however, we will now present some examples of 
vacua which give rise to two strongly coupled hidden sectors. 



IV. EXAMPLES OF G 2 MANIFOLDS 

Having shown that the potential stabilizes all the moduli, it is of interest to construct explicit examples 
of 02-manifolds realizing these vacua. To demonstrate the existence of a G2-holonomy metric on a 



compact 7-manifold is a difficult problem in solving non- linear equations 211 ] . There is no analogue of 



Yau's theorem for Calabi-Yau manifolds which allows an "algebraic" construction. However, Joyce and 



Kovalev have successfully constructed many smooth examples [21[. Furthermore, dualities with heterotic 
and Type IIA string vacua also imply the existence of many singular examples. The vacua of interest to 
us here are those with two or more hidden sector gauge groups These correspond to G2-manifolds which 
have two three dimensional submanifolds Q\ and O2 along which there are orbifold singularities. In order 
to describe such examples we will a) outline an extension of Kovalev's construction to include orbifold 
singularities and b) use duality with the heterotic string. 

Kovalev constructs G2 manifolds which can be described as the total space of a fibration. The fibres 
are four dimensional K3 surfaces, which vary over a three dimensional sphere. Kovalev considers the 
case in which the K3 fibers are generically smooth, but it is reasonably straightforward to also consider 
cases in which the (generic) K3 fiber has orbifold singularities. This gives G2-manifolds which also have 
orbifold singularities along the sphere and give rise to Yang-Mills fields in M theory. For example if 
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the generic fibre has both an SU(4) and an SU (5) singularity, then the G2 manifold will have two such 
singularities, both parameterized by disjoint copies of the sphere. In this case and Nf are equal 
because Q\ and Q2 are in the same homology class, which is precisely the special case that we consider 
both above and below. 

We arrive at a very similar picture by considering the M theory dual of the heterotic string on a 
Calabi-Yau manifold at large complex structure. In this limit, the Calabi-Yau is T 3 fibered and the 
M theory dual is -KT3-fibered, again over a three-sphere (or a discrete quotient thereof). Then, if the 
hidden sector E% is broken by the background gauge field to, say, SU (5) x SU (2) the _RT3-fibers of the 
Cr2-nianifold generically have SU(5) and SU(2) singularities, again with = Nf. More generally, in 
K3 fibered examples, the homology class of Qi could be k times that of Q2 and in this case = kNf. 
As a particularly interesting example, the M theory dual of the heterotic vacua described in [22| include 
a G2 manifold whose singularities are such that they give rise to an observable sector with precisely the 
matter content of the MSSM whilst the hidden sector has gauge group G = Eg. 

Finally, we also note that Joyce's examples typically can have several sets of orbifold singularities 
which often fall into the special class [2l|]. We now go on to describe the vacua in which supersymmetry 
is spontaneously broken. 



V. VACUA WITH SPONTANEOUSLY BROKEN SUPERSYMMETRY 

The potential Q also possesses vacua in which supersymmetry is spontaneously broken. Again these 
are isolated, so the moduli are all fixed. These all turn out to have negative cosmological constant. We 
will see in section IVII that adding matter in the hidden sector can give a potential with de Sitter vacua. 

Since the scalar potential ([9]) is extremely complicated, finding solutions is quite a non-trivial task. 
As for the supersymmetric solution, it is possible to simplify the system of N transcendental equations 
obtained. However, unlike the supersymmetric solution, we have only been able to do this so far for the 
special case as in IIIIBI Therefore, for simplicity we analyze the special case in detail. As we described 
above, there are examples of vacua which fall into this special class. Moreover, as explained previously, 
we expect that typically vacua not in the special class are beyond the supergravity approximation. 
By extremizing (|14p with respect to we obtain the following system of equations 

N 

2ul{bxa - b 2 ) 2 - ffc[2(&ia - 6 2 )(6?a - b\) ^aji/f + 3(a - 1)0? a - b\)v ■ a 

i=l 

N 

+3(6i« — b2) 2 v ■ a + 3(a — l)(b\a — 62)] — 3[(6ia — 6 2 ) 2 a-ivf 

i=l 

+3(a - l)(&ia - b 2 )v- a + 3(a - l) 2 ] =0, (35) 
where we have again introduced an auxiliary variable a defined by 

4?- = I e -(k-&2)*3. (36) 
A\ a 

similar to that in section IIIIBI The definition (|36|) together with the system of polynomial equations 
(f35|) can be regarded as a coupled system of equations for a and vy.. We introduce the following notation: 

x = (a — 1) , y = (ha - b 2 ) , z = (b\a — b 2 ), w = ^ . (37) 
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In this notation, from (|35p (divided by x ) we obtain the following system of coupled equations 

N \ / 2 N \ 

J^ai^ 2 + 3-(u; + l)^-a + 3 -z/ fe - 3 ^ ^ a^ 2 + 3-z7 • a + 3 =0. (38) 




x /x r * — * x 

/ \ i=i / 

It is convenient to recast this system of N cubic equations into a system of N quadratic equations plus a 
constraint. Namely, by introducing a new variable T as 

2 N 

4r = 2^wVa;ZA 2 + 3-(«; + lW-a + 3, (39) 
x l x 
i=i 

where the factor of four has been introduced for future convenience, the system in ([38]) can be expressed 
as 

N 



2 / 2 JV \ 

2 ^" 4T f^" 3 (^^ ail/i2+3 f Z7 ' (?+3 ) (40) 

An important property of the system (|40p is that all of its equations are the same independent of the 
index k. However, since the combination in the round brackets in (|40p is not a constant with respect 
to V this system of quadratic equations does not decouple. Nevertheless, because both the first and the 
second monomials in (|40p with respect to v k are independent of v, the standard solution of a quadratic 
equation dictates that the solutions for i/j. of (|40p have the form 

x 



u k = - (T + m k H) , with m k = ±1 , k = 1, N , (41) 

y 

where we introduced another variable H and pulled out the factor of x/y for future convenience. 

We have now reduced the task of determining v k for each k = 1,N to finding only two quantities - 
T and H. By substituting (|4"Tj) into equations (|39H4Up and using ([5]), we obtain a system of two coupled 
quadratic equations 

^ (Tj + 2AT A H A + H\) + 7(w + 1) (T A + AH A ) + 3 - 4T A = (42) 
9 (Tj + 2AT A ^ + H\) - AH A (H A + AT a) + 21 (T A + AH a) + 9 = 0, 
where parameter A defined by 

3 

A = -m-a, (43) 

is now labelling each solution. Note that by factoring out x/y in (|4ip . the system obtained in (|42p is 
independent of either x or y. However, it does couple to the constraint ([36]) via w. In subsection IV CI we 
will see that there exists a natural limit when the system (|42p completely decouples from the constraint 
([36|) . Since both Ta and both depend on the parameter A, the solution in ([4T]) is now written as 

v£ = -(T A + m k H A ) • (44) 

y 

Since k = X,N and m k = ±1, vector m represents one of 2^ possible combinations. Thus, parameter 
fa ■ a can take on 2 N possible rational values within the range: 



(45) 
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so that parameter A defined in (|43p labelling each solution can take on 2 N rational values in the range: 

-1<A<1. (46) 
For example, when N = 2, there are four possible combinations for m = (mi, 7712), namely 

(mi,m 2 ) = {(-l,-l), (1,-1), (-1,1), (1,1)}. (47) 
These combinations result in the following four possible values for A: 

A = {-1, -(01 - 02) , y(-oi + 02) , 1} , (48) 

where we used ([2]) for the first and last combinations. 

In general, for an arbitrary value of A, system (|42p has four solutions. However, with the exception of 
the case when A = 1, out of the four solutions only two are actually real, as we will see later in subsection 
IV CI The way to find those solutions is the following: 

Having found analytically in terms of a and the other constants, we can substitute it into the tran- 
scendental constraint ([36]) to determine a numerically for particular values of {Ai, A2, b\, 62, N^, a^}. 
Again, in general there will be more than one solution for a. We can then substitute those values back into 
the analytical solution for to find the corresponding extrema, having chosen only those a, obtained 
numerically from (|36p . which result in real values of vfi- We thus have 2 N+1 real extrema. However, after 
a closer look at the system of equations (I42D we notice that when A — > —A, equations remain invariant 
if H-a —> —Ha, and T^a — > Ta, thus simply exchanging the solutions vA + ^ where = 1 with i/^j. _x 
where rrik = — 1, i.e. 

< +) = <_) , (49) 

which implies that the scalar potential (I14j) in general has a total number of 2 N real independent ex- 
trema. However, as we will see later in section IV C\ many of those vacua will be incompatible with the 
super gravity approximation. 

For general values of A, equations (|42p have analytical solutions that are too complicated to be 
presented here. In addition to restricting to the situation with the same gauge kinetic function / in both 
hidden sectors, we now further restrict to special situations where A takes special values, so that the 
expressions are simple. However, it is important to understand that they still capture the main features 
of the general solution. In the following, we provide explicit solutions (in the restricted situation as 
mentioned above) for M theory compactifcations on G2 manifolds with one and two moduli respectively. 
In subsection IV CI we will generalize our results to the case with many moduli and give a complete 
classification of all possible solutions. We will then consider the limit when the volume of the associative 
cycle Vol(Q) = V • a is large and obtain explicit analytic solutions for the moduli. 

A. One modulus case 

The first, and the simplest case is to consider a manifold with only one modulus, i.e. iV = 1, a = I. 
In this case, A = ±1. From the previous discussion we only need to consider the case A = 1. It turns 
out that this is a special case for which the system (j42j) degenerates to yield three solutions instead of 
four. All three are real, however, only two of them result in positive values of the modulus: 

tf» = ~. flf'-l «») 
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and 



3 3645 

T\ 2) = — — (_i34i9 + + 15288k; - 5488kj 2 

1 28 (243 - Ulw + 196kj 2 ) v w 

-329\/729 - 1701k; + 1323kj 2 - 343k; 3 (51) 
135 



+ V729 - 1701u; + 1323kj 2 - 343kj 3 

w 

+196^^729 - 1701k; + 1323kj 2 - 343kj 3 ) 

H m = JL C_27 + 28kj - y/729 - 1701k; + 1323k; 2 - 343k; 3 ) , 
28k; V / 

which give the following two values for the modulus 

s (2) = M 2) + h[ 2) ) = -— ^ + ^ 9 ~ 7w 



iVi y \ 1 1 J Ny \ w 

In addition, each solution in (|52|) is a function of the auxiliary variable a defined in (|36|) . By substituting 
(|52p into (|36p we obtain two equations for and 

Ai " ad) ' A 1 ~ «( 2 ) 6 " {M) 

The transcendental equations (|53p can only be solved numerically. Here we will choose the following 
values for this simple toy model 

2-7T 2ir 

Ai = 0.12, A 2 = 2, 6i = — , b 2 = —, N 1 = l. (54) 

O I 

By solving (|53p numerically and keeping only those solutions that result in real positive values for the 
modulus s in (f52l) we get 

s (1) = 26.101 , s (2) = 27.185 . (55) 

with 

a (1) = 1.122, = 1.267. (56) 

In figure © we see that the two solutions in (|55p correspond to an AdS minimum and a de Sitter 
maximum. In fact, the AdS minimum at is super symmetric. The general solution for given in 
(|52p can also be obtained by methods of section IIIIl imposing the SUSY condition on the corresponding 
F-term by setting it to zero, while introducing the same auxiliary constraint as in (|36l) . 



B. Two moduli case 



While the previous example with one modulus is interesting, it does not capture some very important 
properties of the vacua which arise when two or more moduli are considered. In particular, in this subsec- 
tion we will see that the supersymmetric AdS minimum, obtained in the one-dimensional case, actually 
turns into a saddle point whereas the stable minima are AdS with spontaneously broken super symmetry. 
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FIG. 6: Potential multiplied by 10 32 plotted as a function of one modulus s. For our particular choice of constants in 
(|54l) . the modulus is stabilized at the supersymmetric AdS minimum s' 1 ' = 26.101. The maximum is de Sitter, given by 
s (2) = 27.185. 

Let us now consider a particularly simple example with two moduli. Here we will choose both moduli to 
appear on an equal footing in the Kahler potential ((T|) by choosing 

7 7 

ai = - , a 2 = - . (57) 
We now have four possible combinations for m = (7711,777-2): 

(1,1), (1,-1), (-1,1), (-1,-1), (58) 
corresponding to the following possible values of A: 

1, 0, 0, -1, (59) 

where only two of the four actually produce independent solutions. The case when A = 1 has been solved 
in the previous subsection with t[ 1 \ h[ 1} and t[ 2 \ h[ 2) given by (|50H5ip with the moduli taking on the 
following values for the supersymmetric AdS extremum 

(i) a\x ( 3\ 7x 



s 



(1) _ Q-2X ( 3\ 7x 



(60) 



2 



N 2 y V 2 J AN 2 y ' 



and the de Sitter extremum 



w 

2 N 2 y\7w y 7 2N 2 y \ w 

As mentioned earlier, the supersymmetric solution can also be obtained by the methods of section IIII1 
Now, we also have a new case when A = 0. The corresponding two real solutions for To and ^0 are 

r ° (1)= iiL^ 5 " 63 ^"^' (62) 

/i) 1 / 585 18225 3915 1215 ~ 225 ~ 

Hk ' = — -=\ tt H 1 ttD D , 

4^/5 V 8 39277J 2 28w 392w 2 56w ' 
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and 



T ° (2) = iiL (15 " 63u, "- D) ' (fi;!) 



„(2) 1 / 585 18225 3915 1215 „ 225 „ 

Hk = =W tt H 1 T D D 

4^/5 V 8 392w 2 28w 392w 2 56w 



where we denned 



D = V225 - 770w + 833-u; 2 . (64) 
The moduli are then extremized at the values given by 

•P-^( s f )+| tf ) )' < 65 > 

(3) _ a 2 x / (i) (i) 

2 ~n^\° ~ ° 



and 



'S 4, =ii( T « (2,+^^ ^ 2, )• < 66 > 



a 2 z / T (2) „( 2 )\ 



To completely determine the extrema we again need to substitute the solutions given above into the 
constraint equation (|36p and choose a particular set of values for A\, A2, b\, and 62 to find numerical 
solutions that result in real positive values for the moduli s% and S2- Here we again use the same values 
as we chose in the previous case given by 

2-7T 2ir 

Ax = 0.12, A 2 = 2, h = —, b 2 = —, Ni = l, N 2 = l. (67) 



For the SUSY extremum we have 



The de Sitter extremum is given by 



= 13.05, = 13.05. 



4 2) = 13.59 , s 2 2) = 13 - 59 • (69) 



The other two extrema are at the values 



4 3) = 2.61 , 4 3) = 23 - 55 . and s i 4) = 23 - 55 > 4 4) = 2 - 61 ( 70 ) 



It is interesting to note that the super symmetric extremum in (|68p is no longer a stable minimum 
but instead, a saddle point. The two symmetrically located stable minima seen in figure ([7]) are non- 
supersymmetric. Thus we have an explicit illustration of a potential where spontaneous breaking of 
supersymmetry can be realized. The stable minima appear symmetrically since both moduli were chosen 
to be on an equal footing in the scalar potential. With a slight deviation where a% 7^ a 2 and/or N\ 7^ N 2 
one of the minima will be deeper that the other. It is important to note that at both minima, the volume 
given by ([2]) is stabilized at the value Vx = 122.28 which is large enough for the supergravity analysis 
presented here to be valid. 
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FIG. 7: Potential multiplied by 10 32 plotted as a function of two moduli si and S2 for the values in (|67[) . The SUSY AdS 
extremum given by (|68|l is a saddle point, located between the non-supersymmetric AdS minima given by J7DJ. 



C. Generalization to many moduli 



In the previous section we demonstrated the existence of stable vacua with broken SUSY for the special 
case with two moduli. Here we will extend the analysis to include cases with an arbitrary number of 
moduli for any value of the parameter A. It was demonstrated in section Ull Bl that the SUSY extremum 
has an approximate analytical solution given by (I29p . Therefore, it would be highly desirable to obtain 
approximate analytical solutions for the other extrema in a similar way. We will start with the observation 
that for the SUSY extremum (|60p obtained for the special case when A = 1, both and h[ 1} given 
by (|50jl are independent of w. On the other hand, if in the leading order parameter a is given by (128|) . 
from the definitions in (|37p it follows that in this case 

y — > , and w — ► — oo . (71) 

Thus, if we consider the system (|42p in the limit when w — ► — oo, we should be able to still obtain the 
SUSY extremum exactly. In addition, one might also expect that the solutions for the vacua with broken 
SUSY may also be located near the loci where y — > 0. With this in mind we will take the limit (|7ip 
which results in the following somewhat simplified system of equations for Ta and Ha'- 

2 (Tj + 2AT A H A + H A ) + 3 (T A + AH a) = (72) 
9 [Tl + 2AT A H A + H\) - 4H A (H A + AT a) + 21 (T A + AH a) + 9 = 0. 

Because system (I72D is completely decoupled from the constraint (|36p and hence the microscopic con- 
stants, we can perform a completely general analysis of the vacua valid for arbitrary values of the mi- 
croscopic constants, at least when the limit (|7ip is a good approximation. It is straightforward to see 
that (|50p is an exact solution to the above system when A = 1. Moreover, unlike the general case when 
w is finite, where the system had three real solutions two of which resulted in positive moduli, system 
(|72p above completely degenerates when A = 1 yielding only one solution corresponding to the SUSY 
extremum. On the other hand, for an arbitrary < A < 1 the system has four solutions. One can check 
that at every point A in the range < A < 1 exactly two out of these four solutions are real. The 
corresponding plots are presented in Figure [H Before we discuss the plots we would like to introduce 
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some new notation: 



L { 2\ = Ti c) + m k Hf , (73) 



where c = 1, 2 corresponding to the two real solutions. In this notation (|44p can be reexpressed as: 

■ft = ^2* ■ ( 74 ) 

The volume of the associative three cycle Q for these vacua is then: 

N N 

7 i c) = Val(Q)S> = = £ *< s% = £ * "ft = ^ 3 • 4 C) • (75) 

t=i i=i y 



For future convenience we will also introduce 

'(c) = 

3 



< = a-4 C) = J W C) +M C) ) • (76) 



Constraint (|36p is then given by: 



which is coupled to 



^J^.-M^, (77) 



(c ) -, /r- (c) 



x aV — 1 7 



y ha i 2 ) -b 2 B {c) 



(78) 



.4 



fc) (c) 

where definitions ([37]) were used to substitute for cc and y. Both fe and By are completely determined 

by the system (I72p . whereas 7^ is determined from (|77ti78p . Then solution (j74H can be conveniently 
expressed as 



(c) _ 7 A - (c) 



(<0 



' •.'!./, ,A L A,k- ( 79 ) 

(c) 

Recall that m k = ±1. Thus the only two possibilities for L A k for any k = 1, N are 

L f ± = T A C) ± i^ c) . (80) 

As we vary parameter ^4 over the range < A < 1 point by point, system (|72|) always has exactly two 
real solutions. In Figure[8]we present plots of L A ,, L A _ and B A , where c = 1, 2 as functions of A. We 
only need to consider the positive range < A < 1 because of the symmetry (j49]). 

What happens to these solutions when A = 1? We already know from the previous discussion that 
the system H72\i obtained in the limit w — > — oo degenerates for ^4 = 1 and one obtains the solution that 
corresponds to the SUSY extremum explicitly. The solutions plotted in Figure [8] were obtained assuming 
A ^ 1 and therefore have an apparent singularity when A = 1. Thus they cannot capture either the 
SUSY or the de Sitter extrema that arise in this special case. To explain what happens to the de Sitter 
extremum we need to examine the exact solution in (|5ip . in the same limit. Indeed, bearing in mind that 
w is negative, from (|5ip we have 

-(2) _^(2) , „(2) 3 f 3 + V9 - 7w 



L% + = T?>+Hr=-- r ^ ) . (81) 
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FIG. 8: Plots of L A+ , L A _ and B A ° , where c = 1, 2, corresponding to the two real solutions of the system l|72[> as functions 
of parameter A in the range < A < 1. Both left and right graphs have L A - long dashed line, L A _ -short dashed line, 
B^ A - solid line. 

Left: Plots of L A+ , L A _ and B^ A corresponding to the first real solution at each A. There is a critical value A = 1/7 
where L A _ = and becomes positive for A > 1/7. 

(2) (2) f2) 

Right: Plots of L A ,, L A _ and B A corresponding to the second real solution at each A. 



Here we see immediately that in the limit w — > — oo for the solution above L A + — ► 0. Therefore we 
conclude that the de Sitter extremum cannot be obtained from (|72p which correlates with the previous 
observation that for A = 1 (|72j) has only one solution - the SUSY extremum. Nevertheless, as we will see 
in the next subsection the real solutions plotted in Figure [8] are a very good approximation to the exact 
numerical solutions corresponding to the AdS vacua with spontaneously broken super symmetry. 

Now we would like to classify which of these AdS vacua have all the moduli stabilized at positive 
values. Indeed if some of the moduli are fixed at negative values we can automatically exclude such 
vacua from further consideration since the supergravity approximation assumes that all the moduli are 
positive. Since the volume Tjf 1 is always positive by definition, from (|79h we see immediately that for all 

(c) (c) (c) 

moduli to be stabilized in the positive range, all three quantities L A , , L A _ and B A must have the same 
sign. In Figure [8] the plots on the right satisfy this requirement for the entire range < A < 1. On the 
other hand, the short-dashed curve corresponding to L A _ on the left plot is negative when < A < 1/7, 

features a zero at A = 1/7 and becomes positive for 1/7 < A < 1. Yet, both L A+ and B A remain 

negative throughout the entire range. Moreover, it is easy to verify that the solution with Ty 7 = —3/4 

and Hylj = —3/4, such that Ly 7 _ = is also an exact solution for the general case (|42p when w is 
finite. Therefore, all solutions compatible with the SUGRA approximation can be classified as follows: 

Given a set of {a^} with i = 1, N, there are 2^ possible values of A, including the negative ones. From 
the symmetry in (|49p . only half of those give independent solutions. This narrows the possibilities to 
2 N ~ l positive combinations that fall in the range < A < 1. For each A in the range < A < 1/7 there 
exist exactly two solutions describing AdS vacua with broken SUSY with all the moduli fixed at positive 
values. 

For each A in the range 1/7 < A < 1 there exists exactly one solution describing an AdS vacuum with 
broken SUSY with all the moduli stabilized in the positive range of values. For A = 1 there are exactly 
two solutions with all the moduli stabilized in the positive range - de Sitter extremum in (I5ip and the 
SUSY AdS extremum in (|50p . These two solutions are always present for any set of {a{\. 
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D. Explicit approximate solutions 



In this section we will complete our analysis of the AdS vacua and obtain explicit analytic solutions for 
the moduli. We will take an approach similar to the one we employed in section Ull Bl when we obtained 
an approximate formula (|29p . Expressing oS c ^ from (|78p we obtain 

There exists a natural limit when the volume of the associative cycle TS is large. Just like in the 
approximate SUSY case in (|28p . the leading order solution to (I82p in this limit is given by 



0i 



( C ) _ b 2 

1 

independent of A and c. Plugging this into f|77p and solving for we have in the leading order 

T (c) 1 , (A 2 b 2 \ 1 PQ (A 2 P\ 

where we again assumed the hidden sector gauge groups to be SU(P) and SU(Q). Notice that this 
approximation automatically results in the limit w — > — oo and therefore, L^,, ^a- anc ^ com Puted 
by solving (|72p and plotted in Figure [8] are consistent with this approximation. Thus, combining (|84p 
with (|79p and f|23|) we have the following approximate analytic solution for the moduli in the leading 
order: 

( C ) ^ 1 fa k \ ( L [ 2] k \ PQ (A 2 P 




M "T^ • (85) 



P-Q \AiQ 

To verify the approximation we can check it for the previously considered special case with two moduli 
when a\ = a 2 = 7/6, i.e. the case when A = 0. By solving (I72p we obtain: 



r(l) 
^0,+ 


_ r(2) 
- ^0,- 


3 

~ 16 


(-9 + 


vTr + 


r(l) 


_ r (2) 


3 

~ 16 


(-9 + 




^0 


_ jR (2) 




-9 + v 


/ 17) . 



-26 + wVTf) , (86) 



-26+ 10V17 



Thus, we have the following two solutions for the moduli for the AdS vacua with broken SUSY: 



and 



s (i) 
s o,i 



,(i) 



.(2) 



1 V^26 + WVT7\ 1 PQ . ( A 2 P\ 0.016 PQ ( A 2 P\ 

1 9-V17 i^^Q ln Kmq) ~ ~ T^q ln toJ ' (87) 



1 ^-26+10^ ^ 1 PQ ^ /A 2 P\ ^ 0.143 PQ X J MP 



9_^17 UttN 2 P-Q \AxQJ N 2 P-Q V^iQ 



V-26 + IOvTA 1 PQ , ( A 2 P\ 0.143 PQ , ( A 2 P 



So ' 1_ l 1+ 9-VT7 ) A*N,P-Q lU \lhQ)~ N 2 P-Q ln y Al Q 



( 2 ) _ ( ^-26 + 10^ ^ 1 PQ (MP\ 0-016 PQ / M 2 P 
''°' 2 "1 9-VT7 J 4vriV 2 P-Q n UiQ/ N l P ~Q \ A iQ 



0.016 


PQ 




P-Q 


0.143 


PQ 




P-Q 


0.143 


PQ 


N 2 


P-Q 


0.016 


PQ 
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The choice of the constants given in f)6T[) results the following values: 

4?i = 4*1 = 2 -62 , 4*2 = 4*1 = 23-64 . (89) 

A quick comparison with the exact values in (|70p obtained numerically leads us to believe that the 
approximate analytical solutions presented here are highly accurate. This is especially true when the 

(c) 

volume of the associative cycle is large. For the particular choice above the approximate value is: 

T (1) = T (2) = 26.16, (90) 

which is indeed fairly large. To complete the picture, we also would like to include the first subleading 
order contributions to the approximate solutions presented here. After a straightforward computation 
we have the following: 

(o P M c) (P-Q\ 2 

and 



Q i n (^p) V Q 




By combining (|92p with (|79p and (|23p it is easy to obtain the corresponding expressions for the moduli 
that include the first subleading order correction: 

( C ) 1 ( *k\ ( L Ak\ PQ , M 2 P\ 4% / a fc \ / p-q , 



VI. VACUA WITH CHARGED MATTER IN THE HIDDEN SECTOR 



Thus far, we have studied in reasonable detail, the vacuum structure in the cases when the hidden 
sector has two strongly coupled gauge groups without any charged matter. It is of interest to study how 
the addition of matter charged under the hidden sector gauge group changes the conclusions. We argue 
that the addition of charged matter can give rise to Minkowski or metastable de Sitter (dS) vacua due to 
additional F-terms for the hidden sector matter fields. Hence, dS vacua are obtained without adding any 
anti-branes which explicitly break super symmetry. This possibility was first studied in j23|. Moreover, 
we explain why it is reasonable to expect that for a given choice of G^-manifold, the dS vacuum obtained 
is unique. 



A. Scalar Potential 



Generically we would expect that a hidden sector gauge theory can possess a fairly rich particle 
spectrum which, like the visible sector, may include chiral matter. For example, an SU{N C ) gauge theory 
apart from the "pure glue" may also include massless quark states Q and Q transforming in iV c and N c 
of SU(N C ). When embedded into M theory the effective superpotential due to gauging condensation for 
such a hidden sector with Nf (Nf < N c ) quark flavors has the following form 24|: 



W = At e Nc ~ N f ^ l=1 1 det(QQ) Nc ~ N f . (94) 
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We can introduce an effective meson field <f> to replace the quark bilinear 

1 /2 

<P = (2QQ) = <p e w , (95) 

and for notational brevity we define 

h = 27T , a = . (96) 

N c -N f ' N c — Nf y ! 

Here we will consider the case when the hidden sector gauge groups are SU(N C ) and SU(Q) with Nf 
flavors of the quarks Q (Q) transforming as iV c (N c ) under SU(N C ) and as singlets under SU(Q). In this 
case, when Nf = 1, the effective nonperturbative superpotential has the following form: 

W = A 1( p a e ibl f + A 2 e ib2 f . (97) 

One serious drawback of considering hidden sector matter is that we cannot explicitly calculate the 
moduli dependence of the matter Kahler potential. Therefore we will have to make some (albeit reason- 
able) assumptions, unlike the cases studied in the previous sections. In what follows we will assume that 
we work in a particular region of the moduli space where the Kahler metric for the matter fields in the 
hidden sector is a very slowly varying function of the moduli, essentially a constant. This assumption 
is based on the fact that the chiral fermions are localized at point-like conical singularities so that the 
bulk moduli Sj should have very little effect on the local physics. In general, a singularity supporting 
a chiral fermion has no local moduli, since there are no flat directions constructed from a single chiral 
matter representation. Our assumption is further justified by the M theory lift of some calculable Type 
IIA matter metrics as described in the appendix. It is an interesting and extremely important problem 
to properly derive the matter Kahler potential in M theory and test our assumptions. 

Thus we will consider the case when the hidden sector chiral fermions have "modular weight zero" 
and assume a canonically normalized Kahler potential. The scalar potential is invariant under Q <-> Q 
and Q = Q along the .D-flat direction. For the sake of simplicity, we will first study the case Nf = 1, but 
later it will be shown that all the results also hold true for Nf > 1. The meson field = (2QQ) 1 / 2 along 
the D-flat direction is such that the corresponding Kahler potential for 4> is canonical. The total Kahler 
potential, i.e. moduli plus matter thus takes the form: 

K = -3 ln(4vr 1 / 3 V x ) + Q j Q + Q ] Q = -3 ln(4vr 1 / 3 V x ) + <f>4> ■ (98) 

The moduli F-terms are then given by 

F k = ie^^NkihA^e'^-^ 1 -^-^ 

+ ^^a e -b 1 N-s+i(b 1 -b 2 )N-t+ia0 + A^N-s^ _ (g9) 

2s k 

In addition, an F-term due to the meson field is also generated 



Fa, = (bne 



-ie+ib 2 N-t 



± + A^a^N-s^-b^N-t+iae + A^N-s 

'o 



(100) 
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The supergravity scalar potential is then given by: 



<t> 2 

407TK S 



N 

X 

8=1 



X 



Y, a Mf + 3(i7- a)(6iAf^ a e- 2bli7 - ff + M|e" 262i7,ff +(&i + fo2Mi^ 2 ^e~ (6l+62)i7 - 5 

8=1 

cos((6i - 62)-^ ■ «#)) + 3(Af 0§ a e - 26li; - s + A\e- 2h ^- s + 2^i^ 2 ^ e -( f ' 1+62 ) !7 - s (101) 



xcos((5i - 6 2 )iV • t + ad)) + ^4 (A 2 4 a f 4 + 1 ) e" 2 ^ + yl 2 ( r< 7 





+2^i A 2 0g ^ + 1 j e-^ +b ^- ff cos((&i - 6 2 )7V • t + o0))] . 

Minimizing this potential with respect to the axions and 9 we obtain the following condition: 

sin((&! - b 2 )N ■ t + aO) = . (102) 
The potential has local minima with respect to the moduli s, when 

cos((6i - b 2 )N ■ t + aO) = -1 . (103) 

In this case (llOip reduces to 

<f> 2 N 
V = -£^[(b 1 A 1 r e- b ^ S -b 2 A 2 e- b ^n 2 Y,aM) 2 

x i=l 

+3(u- a)(Ai#}e- ftl * s - A 2 e- b ^ d )(b 1 A 1 ^e- bl ^ d - b 2 A 2 e~ b ^ d ) (104) 
+3(A 1 ^e- b ^ s - A 2 e- b ^ s ) 2 + -AA^l + fa) e~ b ^ s - A 2 fae~ b ^ d ) 2 \ . 

B. Supersymmetric extrema 

Here we consider a case when the scalar potential f j 1 04 [) possess SUSY extrema and find approximate 
solutions for the moduli and the meson field vevs. Taking into account (|103j) and setting the moduli 
.F-terms (j99|) to zero we obtain 

3 a — 1 

v k = v = --— — , 105) 

2 b\a — b 2 

together with the constraint 

« = ^<t>%e~l ( - bl - b2 > . (106) 
A 2 

At the same time, setting the matter F-term (jlOOp to zero results in the following condition: 

-^ + 1)5-1 = 0. (107) 

n J 
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Expressing d from f)105|) and substituting it into (|107p we obtain the following solution for the meson 
vev at the SUSY extremum: 



2 6 2 + 3/(2^ 



Recall that in our analysis we are considering the case when P = N c — Nf > 0, which implies that 
parameter a defined in (196j) is negative. Thus, since the left hand side of (|108h is positive, for the SUSY 
solution to exist, it is necessary to satisfy 

b 2 >h => P>Q. (109) 

Recall that for the moduli to be positive, the constants have to satisfy certain conditions resulting in two 
possible branches (|30H . Therefore, condition (|109j) implies that the SUSY AdS extremum exists only for 
branch a) in (|30p . In the limit, when v is large, the approximate solution is given by: 

a = £, (no) 

mu 3 PQ (A 2 P 

Si = -r— , with v = — — — — — in 



Ni ' UttP-Q VAQ 

,2 2 



+ 



P-Q P]n 



A 2 P 
AiQ 



where we also assumed that P ~ O(10), such that 0g ~ 1. For the case with two moduli where 
o-i = «2 = 7/6 and the choice 

2-7T 27T 

Ax = 4.1, A 2 = 30, h = —, b 2 = — , N 1 = l, N 2 = 1 , (111) 

the numerical solution for the SUSY extremum obtained by minimizing the scalar potential (I104p gives 

ai«44.5, s 2 ~44.5, O ~ 0.883, (112) 

whereas the approximate analytic solution obtained in (|110p yields 

si«45.0, s 2 ~45.0, 0o~O.882. (113) 

This vacuum is very similar to the SUSY AdS extremum obtained previously for the potential arising 
from the "pure glue" Super Yang-Mills (SYM) hidden sector gauge theory. Thus, we will not discuss it 
any further and instead move to the more interesting case, for which condition (|109[) is not satisfied. 

C. Metastable de Sitter (dS) minima 

Below we will use the same approach and notation we used in section |V1 to describe AdS vacua with 
broken SUSY. Again, for brevity we denote 

x = (a — 1) , y = (bid — b 2 ) , z = (b^a — b 2 ) , w = %~ . (114) 

yZ 
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Extremizing (|104p with respect to the moduli Sj and dividing by x 2 we obtain the following system of 
coupled equations 



2 ^-( 2 |,g^ + 3| ( , + 1) ,. 3+3 + ^(^ + 1 )(|| + 1 ))^ «H5) 
plus the constraint (|106|) . Next, we extremize (|104p with respect to </>o and divide it by 2cJ)qx 2 to obtain: 



V 2 x— -v 9 ^ ad /a — 1 \ l0 \ f aa \ 

|iE^+2r-+4(^- a+ -(^- +2 ) +5+ *«)(*P + 1 ) (116> 

To solve the system of A" cubic equations (|115p . we introduce a quadratic constraint 

Af^2^wj2a i u! + si(w + l)u-a + 3 + h 2 (^ : + l) + 1) , (117) 

x ~i x 2 \n x ) \Kv J 

such that (|115p turns into a system of iV coupled quadratic equations: 

2^vl-Aflv k -?, [ |J Va iI /? + 3? I 7-a + 3+ ^ f ^ + ll )=0. (118) 



Again, the standard solution of a quadratic equation dictates that the solutions for v k of (|118p have the 
form 



i/ A = - (T + m fc #) , with m k = ±1 , k=l,N . (119) 



We have now reduced the task of determining v k for each k = 1, N to finding on/y too quantities - T and 
i7. By substituting (|119p into equations ()1 16H118}) and using (J2j) , we obtain a system of three coupled 
equations 

1 (ri + 2AT A H A + H%)+1 (f A + AflA) + i (| (T A + AH A ) + ^ + 2 ) + 5 + ^ (120) 

x ( - + 1) + * ^ ( (is + ™ + ^ ) ; » (TA + ^ ) = o 

^(^ + 2^^ + J^)+7(«, + l)(f A + ^)+3 + |^ + (|| + l)-4T A = 

9 (ri + 2^T a ^a + Ha) ~ ^ A (h a + AT A ) + 21 (t a + AH A ) + 9 + ^ + 1^ = , 

plus the constraint (|106j) . Note that each solution is again labelled by parameter A so that (|119p becomes 

v£ = % (f A + m k H A ) . (121) 
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Let us consider the case when A = 1. In this case, the solution is given by 



v=%(f x + Ht) = %L lt + . (122) 

y v / y 



and (|120p is reduced to 

l(f 1+ H 1 )\l(f 1 + H 1 ) + 3 1 (^(f 1 + H 1 ) + f(^ + 2) + 4 + 5)(^ + l) (123) 



aabi 7 / / ~ ~ \ 2 3 , .., , . 



3y V v ' 2 



14w ^ 1 + ^0 a + 7(«, + i)(r 1 + J & 1 ) + 3 + ^ f^ + l) (§| + l)-4r 1 = 

2 _ / _ _ \ / _ \ vi . / aa 



9[T! + Hi) -4H 1 ^H 1 + T 1 )+21[T 1 +H 1 j+9 + -4> 2 ^-^ + lJ =0. 
In the notation introduced in (|114f) . the SUSY condition (| 107[) can be written as 

^ + x = 0. (124) 

n 

It is then straightforward to check that in the SUSY case, the system (|123|) yields 

as expected. We will now consider branch b) in (|30p for which (|124p is not satisfied. Moreover, in order 
to obtain analytical solutions for the moduli and the meson vev 0o we will again consider the large three 
cycle volume approximation. Recall that in this case we take y — > and w — > — oo limit to obtain the 
following reduced system of equations when A = 1 for Li + and <^o : 

7/~ \ 2 7~ 3/14- aa fa-1 \ 9 \ ( aa \ . , 

J M + j% + 2 ( 3 + - ("jj- + 2) + + 5) + 1 j (126) 

aa6 x 7 //- \2 3~ \ n 
2 /~ \ 2 ~ 3aab\y ( aa \ 

Note that in (|126p . we have dropped the third equation since for A = 1 we only need to know Li.+ and the 
third equation in (|123l) determines Hi + in terms of We also kept the first subleading term in the 

second equation. Note that the term in the second line of the first equation proportional to ~ 1/y appears 

to blow up as y — > 0. However, from the second equation one can see that the combination ^Li i+ ^ 

is proportional to y which makes the corresponding term finite. By keeping the subleading term in the 

second equation, we can express 



2 



) =-^i, + -=^(^ + l) (127) 



j 2 3^ 9aabiy / aa 
2 1,+ ~ 28xz V0p 

from the second equation to substitute into the first equation to obtain in the leading order 



14 ~ 2 aa ( a — \ \ 1 / aab\ \ \ ( aa 



3 r x V 4>q J zx \ 4>q J J \4>qX 



^+ + 5 + <% + -( — + 2 -» Hr H- + 1 )=°- ( 128 ) 
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Since we are now considering branch b) in (|30p . the second factor in (|128p is automatically non-zero. 
Therefore, the first factor in (|128p must be zero. Thus, after substituting 



3 3aabiy ( aa 
2 + Uxz V^p 

obtained from (|127p . we have the following equation for 



L h+ ~ " + ( 72^ + 1 ) . ( 129 ) 



XZ \ X J X \ 0q I ZX 



2 



^H^te + lU^ ^ + 2)4(^) =0. (130) 



Also, since in the leading order L\ t+ = —3/2, using the definitions in (|114j) we can express a from 
(|122p in the limit when f is large, including the first subleading term 

b 2 3(61-62) n „ n 
61 26| ^ 



By combining (|106p with the leading term in (|13ip and taking into account that ^ ~ 1 we again obtain 



= -tt , with 1/ ps — -In -— . (132) 



Thus, from (jl31f) we have 



P 7(Q-P) 2 , , 

a^7T+ /Ax ■ (133) 



Finally, using (|133p along with the definitions of x, y and z in (|114p in terms of a we can solve for 0q 
from (|130p and assuming that Q — P ~ 0(1), in the limit when P is large we obtain 



, 2 . * /. 7/3 



^-^^-^-^^^yi^v 1 -^;- (134> 

We notice immediately that since </>q is real and positive it is necessary that 

Q-P>2. (135) 

We will show shortly that the extremum we found above corresponds to a metastable minimum. Also, 
for a simple case with two moduli, via an explicit numerical check we have confirmed that if Q — P < 2 
the local minimum is completely destabilized yielding a runaway potential. Also note that for (|134p to 

be accurate, it is not only P which has to be large but also the product P In nas ^° s ^ a y l ar § e to 

keep the subleading terms suppressed. To check the accuracy of the solution we again consider a manifold 
with two moduli where the values of the microscopic constants are: 

a 1 = a 2 = 7/6, P = 20, Q = 23, A x = 27, A 2 = 2, N x = N 2 = 1. (136) 

The exact values obtained numerically are: 

si « 33.470 , s 2 ps 33.470 , <p Q ps 0.810 . (137) 
The approximate equations above yield the following values: 

si pa 33.463 , s 2 « 33.463 , <p Q ps 0.803 . (138) 
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Note the high accuracy of the leading order approximation for the moduli s,. 

It is now straightforward to compute the vacuum energy using the approximate solution obtained 
above. First, we compute 

_ 2P 

K^F.F- - Z\W\ 2 = 4 (A 2 x-) 2 (1 ( Lli+ ) 2 + 7 -L h+ + l) Q ~" (139) 

and 



2P 



^ = (W (* +1 )'(*«)-'-' m 

Using (|129p . (|139p and (|140p we obtain the following expression for the potential at the extremum 
with respect to the moduli s» as a function of 0o 



vn = o- 

64ttV| 



, 2aa \ , 9 f aa\ 



21 ^fMV^ (14 i) 



where the terms linear in y cancelled and the quadratic terms were dropped. A quick look at the structure 
of the potential (|141|) as a function of 0q, where 0g > 0, is enough to see that there is a single extremum 
with respect to </>q which is indeed, a minimum. The polynomial in the square brackets is quadratic 
with respect to 0g. Moreover, the coefficient of the cj>Q monomial is equal to unity and therefore is 
always positive. This implies that for the minimum of such a biquadratic polynomial to be positive, it is 
necessary for the corresponding discriminant to be negative, which results in the following condition: 

3-4— <0. (142) 

x 

Again, since in the leading order £14- = —3/2, using the definitions in (|132|) . we can express a from (|122[) 
in terms of v to get 

a a P 3PQ . , 

- = = 1 . (143) 

x 5-1 P-Q 4vri/(P - Q) y ' 

We then substitute v from (I132p into ()143p and use it together with a = —2/P we obtain from (|142p the 
following condition 

8 28 



<0. (144) 



The above equation is the leading order requirement for the energy density at the minimum to be positive. 
It is also clear that the minimum is metastable, as in the decompactification limit (Vx — ► 00), the scalar 
potential vanishes from above, leading to an absolute Minkowski minimum. Figure [9] shows the scalar 
potential for a manifold with two moduli along the slice s± = S2 with the meson field 4> equal to its value 
at the minimum of the potential (|134|) . The microscopic constants are the same as in f|136|) . 



D. The uniqueness of the dS vacuum 

In the previous subsection we found a particular solution of the system in (I120p corresponding to 
A = 1. Here we would like to investigate if solutions for < A < 1 are possible when the vacuum for 
A = 1 is de Sitter. Just like for the pure Super Yang-Mills (SYM) case, we can recast (11211) as 
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V/ (m 3 / 2 2 m p 2 ) 
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Si=S 2 



FIG. 9: Potential in units of m 



3/2 7 



along the slice si = S2 for a manifold with two moduli with the meson field equal to 



its value at the minimum of the potential (|134[) . The microscopic constants are as in (|136[) . Although hard to see from the 
graph, the value of the potential at the minimum (i.e. the cosmological constant) is 0.194my 2 mp. 



where the volume of the associative three cycle Q is again 

T A = Vol(Q) A = a-P A =*B A , (146) 

y 

and we have introduced 

B A = a-L A = 7 -(f A + AHj^j . (147) 

Just like we did in equation (|82p for the pure SYM case, we can also express a A as 

b 2 T A -B A 

aA = T^P W ■ (148) 

biT A - B A 

If we again consider the large associative cycle volume limit and take y — » and w — * — oo, the second 
and third equations in (|120j) in the leading order reduce to 

2 (f% + 2Af A H A + H 2 A )+3 (f A + AH A ) = (149) 
9 (f 2 A + 2Af A H A + h£) - AH A (h a + AT A ^j + 21 (f A + AH A ^j + 9 + (j^z + lj = . 

Note that the only difference between (|72p and (I149P is the presence of the term 

9 2 / a & 



which couples the system (|149p to the first equation in (|120p which determines 4>q. Instead of solving the 
full system to determine T A , H A and 4>o and analyzing the solutions we choose a quicker strategy for our 
further analysis. Namely, we can solve the system of two equations in (|149p and regard 5 as a continuous 
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deformation parameter. One may object to this proposition because a and therefore x = a — 1 are not 
independent of parameter A. However, in the limit when Ta is large, we notice from (|148p that in the 
leading order, ola is indeed independent of A. 

Recall that in the pure SYM case the system (I72p corresponding to the case when 5 = has two real 
solutions for all < A < 1. Thus, one may expect that as we continuously dial 5, the system may still 
yield real solutions for A < 1. Let us first determine the range of possible values of parameter 5. A quick 
calculation yields that the combination in (|150p is the smallest with respect to (j>o when 0q = ^S-. In this 
case 

aa 

5 = 9—. 151) 

x 

Now, recall from the previous subsection that for the solution corresponding to A = 1 to have a positive 
vacuum energy, condition (|142p must hold. Since a and x are independent of A in the leading order, 
condition (|142p implies that 

27 

5 > (152) 
4 

Again, since the volume Ta is always positive, from (|145p we see that for all moduli to be stabilized in 
the positive range, all three quantities La, + , La,- and Ba must have the same sign. For 5 = 27/4, 
the system (|149p has two real solutions when 0.877781 < A < 1. However, from the left plot in Figure 




FIG. 10: Plots of L A L A _ and B A , where c = 1,2, corresponding to the two real solutions of the system (|149p as 
functions of parameter A. L A , - long dashed curve, L A _ -short dashed curve, B { A ] - solid curve. Black color: 1 '- 

and B A corresponding to the first real solution. Red color: L A + , L A _ and B A corresponding to the second real solution. 
Left plot: when 8 = 27/4 the real solutions exist in the range 0.877781 < A < 1. Right plot: when S = 8 the real solutions 
exist st in the range 0.915342 < A < 1. 



[10] corresponding to the minimum value 5 = 27/4 we see that neither of the two solutions satisfy the 
above requirement since both short-dashed curves corresponding to La - for the two solutions are always 
positive for the entire range 0.877781 < A < 1, whereas both La,+ and Ba remain negative. Therefore 
for 5 = 27/4 and A < 1 there are no solutions for which all the moduli are stabilized at positive values. 
Moreover, as parameter 5 is further increased, the range of possible values of A for which the system 
has two real solutions gets smaller and more importantly, the values of La,- remain positive and only 
increase while both La. + and Ba remain negative, which can be seen from the right plot in Figure [TUl 
where 5 = 8. This trend continues as we increase 5. 

Thus, we can make the following general claim: If the solution for A = 1 has a positive vacuum energy, 
condition (I152p must hold. When this condition is satisfied the system (I149P has no solutions in the range 
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< A < 1 for which all the moduli are stabilized at positive values. Therefore, if the vacuum found for 
A = 1 is de Sitter it is the only possible vacuum where all the moduli are stabilized at positive values. 
Although the above analysis was done in the limit when T4 is large, we have run a number of explicit 
numerical checks for a manifold with two moduli and various values of the constants confirming the above 
claim. In addition, although we have not proved it, it seems plausible from many numerical checks we 
carried out that it is also not possible to have a metastable dS minimum for values of A different from 
unity, even if the dS condition on the A = 1 vacuum is not imposed. 

Finally, it should be noted that the situation with a "unique" dS vacuum is in sharp contrast to that 
when one obtains anti-de Sitter vacua, where there are between 2 Ar_1 and 2^ solutions for N moduli 
depending on the value of A (see section |V|). Let us explain this in a bit more detail. Since the dS 
solution found for A = 1 is located right in the vicinity of the "would be AdS SUSY extremum" 3 where 
the moduli F-terms are nearly zero, it is the large contribution from the matter F-term (I140P which 
cancels the — 3|IU| 2 term in the scalar potential resulting in a positive vacuum energy. Recall that in the 
leading order all the AdS vacua with the moduli vevs s^f • are located within the hyperplane 4 

g = hl^Q 1 " (£5) = constant ■ < 153) 

The matter F-term contribution to the scalar potential K^F^F^ evaluated at the same but arbitrary 
4>q is therefore also constant along the hyperplane (|153p . Thus, while the matter F-term contribution 
stays constant, as we move along the hyperplane (1153j) away from the dS minimum, where the moduli 
F-terms are the smallest, the moduli F-term contributions can only get larger so that the scalar potential 
becomes even more positive. This implies that the AdS minima with broken SUSY found in Section |V] 
completely disappear, as the AdS SUSY extremum becomes a dS minimum. 

VII. RELEVANT SCALES 

We have demonstrated above that in fluxless M theory vacua, strong gauge dynamics can generate a 
potential which stabilizes all the moduli. Since the entire potential is generated by this dynamics, and 
the strong coupling scale is below the Planck scale, we also have a hierarchy of scales. In this section we 
calculate some of the basic scales in detail. In particular, the gravitino mass, which typically controls the 
scale of supersymmetry breaking is calculated. By uniformly scanning over the constants (N,P,Q,Ak) 
with Ni order one, we demonstrate in lVII Cl that a reasonable fraction of choices of constants have a TeV 
scale gravitino mass. We do not know if the space of GVmanifolds uniformly scans the (P,Q,Ak,Ni) 
or not, and more importantly, the scale of variation of the A^'s in the space of manifolds is not clear. 
The variation of the A^s is the most important issue here, since one can certainly vary P and Q over 
an order of magnitude. We begin with a discussion of the basic scales in the problem. We will begin 
with the AdS vacua, then go on to discuss the de Sitter case. In particular, in the dS case, requiring a 
small vacuum energy seems to lead to superpartners at around the TeV scale. It will also be shown that 
including more than one flavor of quarks in the hidden sector or including matter in both hidden sectors 
does not change this result. The section will end with an estimation of the height of the potential barrier 
in these vacua. 



3 This can be seen by comparing the leading order expression for the moduli vevs in the dS case (1 132 p with the corresponding 
formula for the SUSY AdS extremum (|3Tj) . 

4 c=l,2 labels the two real solutions of the system ([72}. 
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A. Scales: AdS Vacua 

As an example, we consider one of the non-SUSY minima in our toy model given by (|70p and compute 
some of the quantities relevant for phenomenology. Namely, the vacuum energy 
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-(5.1 x 10 iU GeV) 



the gravitino mass 



the 11-dimensional Planck scale 
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3/2 



m p e K/2 \W\ 
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7rm ri 
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2.081 TeV . 



3.9 x 10 17 GeV, 
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the scale of gaugino condensation in the hidden sectors 



A (i) 

9 

A (2) 
9 



m p e-^' N ' sl « 2.6 x 10 15 GeV 



9.7 x 10 14 GeV 



(154) 
(155) 

(156) 



(157) 
(158) 



where m p = (8irGx) = 2.43 x 10 18 GeV is the reduced four-dimensional Planck mass. 
From (|155p and (|156p . we see that it is possible to have a TeV scale gravitino mass together with 
Mn > M unif (2 x 10 16 GeV). This feature survives in more general cases as well, implying that standard 
gauge unification is compatible with low scale SUSY in these vacua. 



B. Gravitino mass 



By definition, the gravitino mass is given by: 

™>3/2 = 



m 



n e K / 2 \W\. 



(159) 



For the particular M theory vacua with Kahler potential given by (pQ) and the non-perturbative super- 
potential as in © with SU(P) and SU (Q) hidden sector gauge groups we have: 
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iie -fw_ i2e -|w 



(160) 



where the relative minus sign inside the superpotential is due to the axions. Before we get to the gravitino 
mass we first compute the volume of the compactified manifold Vx for the AdS vacua with broken SUSY. 
By plugging the approximate leading order solution for the moduli (|85p into the definition ([2]) of Vx we 
obtain: 
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(161) 



Recalling the definition (f75l) of Tjf 1 and using (j84l) together with (116ip to plug into (1160p the gravitino 
mass for these vacua in the leading order approximation is given by: 



(^3/2)? 



V2tt 3 AoP 



P-Q 



PQ 



PQ , ( A 2 P 
-in 



P-Q V^iQ 



A 2 P 
A X Q 



N 

n 

i=l 



3aj 

oL (c) 



(162) 



37 



For the special case with two moduli when a\ = a 2 = 7/6, considered in the previous sections we obtain 
the following: 



{m 3/2 ) { o ,2) = m p 2 1 /V (j + Vrry (N!N 2 )* A 2 P 

P-Q 



P-Q 



PQ 



A 2 P\~^Q ( PQ , A 2 P\ a 
m- 



A l Q) \P-Q A t Q 



m p 2.97 x 10 d (A^i N 2 )* A 2 P 



a 2 p\-^q /pq A 2 pyl 

MQ) {p^Mq) (163) 



PQ 

For the choice of constants as in (fBTl) the leading order approximation (1163|) yields: 

(^3/ 2 )i 1,2) = 2061GeV , (164) 
whereas the exact value computed numerically for the same choice of constants is: 

m 3/2 = 2081GeV . (165) 
Again, we see a good agreement between the leading order approximation and the exact values. 

C. Scanning the Gravitino mass 

In previous sections we found explicit solutions describing vacua with spontaneously broken super- 
symmetry. Moreover, we also demonstrated that for a particular set of the constants these solutions 
can result in m^/ 2 ~ O(l) TeV. It would be extremely interesting and worthwhile to estimate (even 
roughly) the fraction of all possible solutions which exhibit spontaneously broken SUSY at the scales of 
O(l)-O(10) TeV. We would first like to do this for generic AdS/dS vacua with a large magnitude of the 
cosmological constant (~ m^^mfy. The analysis for the AdS vacua is given below but as we will see, the 
results obtained for the fraction of vacua are quite similar for the dS case as well. In subsection IVII Dl 
we impose the requirement of a small cosmological constant as a constraint and try to understand its 
repercussions for the gravitino mass. 

We do not yet know the range that the constants (N, P, Q, A\,A 2 ) take in the space of all G 2 manifolds. 
Nevertheless, we do have a rough idea about some of them. For example, we expect that the quantity 
given by the ratio 

(166> 

which appears in several equations, does deviate from unity. One reason for this may be due the threshold 
corrections [191 ] which in turn depend on the properties of a particular GVholonomy manifold. For 
concreteness, we take an upper bound p < 10. Also, based on the duality with the Heterotic String we 
can get some idea on the possible range of integers P and Q corresponding to the dual coxeter numbers 
of the hidden sector gauge groups. Namely, since for both SO(32) and Eg gauge groups appearing in the 
Heterotic String theories the dual coxeter numbers are h v = 30, we can tentatively assume that both P 
and Q can be at least as large as 30. Of course, we do not rule out any values higher than 30 but in this 
section we will assume an upper bound P,Q < 30. 

We now turn our attention to equation (|162p which will be used to estimate the gravitino mass scale. 
It is clear from the structure of the formula that m 3 / 2 is extremely sensitive to P, Q as well as the ratio 
p, given by (|166p . On the other hand it is less sensitive to the other constants appearing in the equation 
such as Ni, di and the ratios /L^\- This is because the powers 3aj/2 for each term under the product 
get much less than one as the number of moduli increases because of the constraint on in ([2]). This will 
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smooth any differences between the contributions coming from the individual factors inside the product. 
Since for < A < 1 (A is defined in (|44p ). the ratios _B^/L^ vary only in the range O(l)-O(10), for our 
purposes it will be sufficient to simply consider (|162|) for the case when A = 1 corresponding to the SUSY 

extremum so that B± /L^) = 7/3 for all i. This is certainly good enough for the order of magnitude 
estimates we are interested in. It also seems reasonable to assume that the integers Ni are all of O(l). 
Yet, even if some Ni are unnaturally large, their individual contributions are generically washed out since 
they are raised to the powers that are much less than one. Thus, for simplicity we will take N = 1 for 



all i = 1,N. Finally, from field theory computations 25], A<i = Q (in a particular RG scheme) up to 
threshold corrections. We therefore take A2 ~ Q for simplicity, allowing A± to vary. 
Thus, the gravitino mass in our analysis is given by 

m 3/2 ~ V2^PQ (^757^) ' M~* fl ' • ( 167 ) 



Finally, with regard to the constants a% which are a subject to the constraint 

N 

cu = 

3 



N 7 



we will narrow our analysis to two opposite cases. For the first case we make the following choice 

1 



1) oi = 2, and ai = ^—- — , for i = 2,N, (169) 

such that one modulus is generically large and all the other moduli are much smaller. This is a highly 
anisotropic GVmanifold. The second case is 

7 



2) a i = — , for all i = l,N, (170) 

with all the moduli being on an equal footing. Therefore, by considering these opposite cases we expect 
that most other possible sets of a, will give similar results that are somewhere in between. For each set 
of a% above, equation (|167h gives 

i) »Si~ ?! ir i «»(^) , w-»w-*(»-i)», <m) 



2) m< 2 / 2 ~ ^PQ * [M _l (Pf 1 ^ . (172) 

For a typical compactification we expect ~ 0(100), therefore the variation of m.3/2 due to an O(l) 
change in the number of moduli for the first case is O(l) whereas in the second case it can be as large as 
O(10). Thus, if we choose N = 100, we expect that our order of magnitude analysis will be fairly robust 
for case 1). For case 2), however, we will perform the same analysis for N = 100 and N = 50 to see how 
different the results will be. Before we proceed further we need to impose a restriction on the possible 
solutions to remain within the SUGRA framework. Using (|16ip . condition that Vx must remain greater 
than one for the two cases under consideration translates into the following two conditions: 



(173) 
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Inp > 1 



(174) 



2vr N P - Q 

Then, as long as conditions (|173M174|) hold, volume of the associative cycle T is greater than one is 
satisfied automatically - a necessary condition for the validity of supergravity. This is obvious from 
comparing the right hand side of (|84|) with each condition above. From (| 173H174j) we can find a critical 
value of p = p cr it for both cases at which Vx = 1: 
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By substituting (|1 75111 76]) into (| 171111721) we can find the corresponding upper limits on 711,3/2 as functions 
of P and Q, below which our solutions are going to be consistent with the SUGRA approximation. 

In Figure [TT1 we present plots of log 10 (ffl.3/2) for both cases as a function of P in the range where 
Pcrit < P < 10 for different values of P — Q. On all the plots the light grey area represents possible 
values of log 10 (m3/ 2 ) consistent with the supergravity framework. For the sake of completeness we have 
also included the formal plot of log 10 (m 3 / 2 ) corresponding to p = 1.01 represented by the dashed curve. 
[From the plots it is clear that as the difference P — Q is increased from 1 to 3 - top and from 1 to 
2 - bottom, both the light grey area representing all possible values of log 10 ("i3/2) consistent with the 
SUGRA approximation and the dark area corresponding to — 16 < log 10 (?H3/2) < — 14 get significantly 
smaller. If we further increase P — Q, the light grey region shrinks even more for case 1), and does not 
exist for case 2), while the dark region completely disappears in both cases. Therefore, for case 2) the 
plots on the bottom of Figure [IT] are the only possibilities where solutions for P < 30 and TV" = 100 
consistent with the SUGRA approximation are possible, implying un upper bound (P — Q) m ax = 2. It 
turns out that for case 1) the upper bound on (P — Q) where such solutions are possible is much higher 
{P Qjmax — 23. 

Assuming that all values of the constants such as P, Q and p are equally likely to occur in the ranges 
chosen above we can perform a crude estimate of the number of solutions with — 16 < log 10 (?H3/2) < — 14 
relative to the total number of possible solutions consistent with the SUGRA approximation. In doing 
so we will use the following approach. For each value of (P — Q) in the range 1 < (P — Q) < (P — Q) ma x 
we compute the area of the grey region for each plot and then add all of them to find the total volume 
corresponding to all possible values of log 10 (?H3/2) consistent with the supergravity approximation. 
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(P—Q)max „3Q 
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(177) 



Likewise, we add all the dark areas for each plot to find the volume corresponding to the region where 
-16 < log 10 (m 3/2 ) < -14 

{P-Q) max p 

n = Yl / dPlog 10 (m 3/2 )| (178) 

(P-Q)=l \ -16<log 10 (m 3/2 )<-14J > 

where (P — Q)* max is un upper bound before the dark region completely disappears. From the previous 
discussion, for case 1): (P — Q)* max = 3; for case 2): (P — Q)* max = 2 . Then, the fraction of the volume 
where 240 GeV < 7113/2 ^ 24TeV is given by the ratio 
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tot 



(179) 
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FIG. 11: log 10 (7713/2) as a function of P for case 1) - top plots and case 2) - bottom plots. The light grey area represents 
possible values of log 10 (m 3/ / 2 ) in the range where p cr u < p < 10 consistent with the SUGRA approximation. The dark area 
indicates the region of interest where — 16 < log 10 (m 3/ / 2 ) < — 14 such that 240 GeV < m 3 / 2 < 24TeV. The dark solid curve 
corresponds to log 10 (m3/ 2 ) when p — p cr u- The lower boundary of the light grey area represents the log 10 (rri3/ 2 ) curve when 
p — 10. The dashed curve corresponds to p — 1.01. Top left: Case 1) when P — Q = 1. Top right: Case 1) when P — Q — 3. 
Bottom left: Case 2) when P — Q = 1. Bottom right: Case 2) when P — Q — 2. 

Numerical computations yield the following values for the two cases when N = 100: 

1) Ai = 3.5%, (180) 

2) A 2 = 13.6%. (181) 

Because of the significant difference in the dependence of pa-it on the number of moduli N in (|175p versus 
(|176h . the number of solutions consistent with the SUGRA approximation is cut down dramatically in 
case 2) compared to case 1). This also occurs because of the different dependence of m 3 / 2 in (|17ip and 
(fl~72l) on N. Namely, for N ~ O(100), the values of m 3/2 for case 2) in (fT72l) are ~ O(10 6 ) greater than 
those for case 1). Furthermore, for the same reasons, it turns out that if we keep increasing the number of 
moduli N, for case 2) there is un upper bound N < 157 for the solutions with 240 GeV < m 3/ / 2 < 24TeV 
compatible with the SUGRA approximation to exist at all. Of course, this upper bound can be higher if 
we allow P to be greater than 30. 

By performing the same analysis for N = 50 we get the following estimates 



1) Ai = 3.4% , 



(182) 
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2) A 2 = 10.7% . (183) 

As expected, decreasing the number of moduli by a half has produced little effect on Ai while decreasing 
A2 by a few present. These numbers coming from our somewhat crude analysis already demonstrate 
that a comparatively large fraction of vacua in M theory generate the desired hierarchy between the 
Planck and the electroweak scale physics. Also, one can easily check that all the solutions consistent 
with the SUGRA framework for which 240 GeV < 7713/2 < 24TeV for any number of moduli N satisfy 
the following bound on the eleven dimensional scale 

3.6 x 10 16 GeV < m u < 4.3 x 10 18 GeV , (184) 

which makes them compatible with the standard unification at Mqut ~ 2 x 10 16 GeV. This is also a 
nice feature. Of course, apart from determining the upper and lower bounds on the constants, it would 
be desirable to know their distribution for all possible manifolds of G 2 holonomy. In this case instead of 
using the flat statistical measure as we did here, each solution would be assigned a certain weight making 
the sampling analysis more accurate. However, this is an extremely challenging task which goes beyond 
the scope of this work. 

The simple analysis presented in this section clearly points to a very restrictive nature of the solutions. 
Namely, the requirement of consistency with the supergravity regime results in very strict bounds on the 
properties of the compactification manifold. Further requirements coming from the SUSY breaking scale 
to be in the range required for supersymmetry to solve the hierarchy problem narrows down the class of 
possible G 2 holonomy manifolds even more. It would be extremely interesting to know to what extent 
these results extend into the small volume, "stringy" regime, about which we have nothing to say here. 



Results for dS Vacua 



We will now show that the results obtained in the previous subsection also hold true for dS vacua 
with Kahler potential given by (I98p and the non-perturbative superpotential as in (|97h with SU(N C ) and 
SU(Q) hidden sector gauge groups. For this case, we have : 
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where the relative minus sign inside the superpotential is due to the axions and P = N c — 1. Before we 
get to the gravitino mass we first compute the volume of the compactified manifold Vx for the metastable 
dS vacuum with broken SUSY. By substituting the approximate leading order solution for the moduli 
(|132p into the definition ([2]) of Vx we obtain: 
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Recalling the definition of Im(/) in terms of v and using the solution for v (eqn. (|132p ) together with 
(|186f) . the gravitino mass for the dS vacuum in the leading order approximation is given by: 
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(187) 



where 4>q is given by fll34H . Since ^q^ 1 ^ ~ 1 from section PVTl and A 2 ~ Q, we see that the expression for 
the gravitino mass for dS vacua is almost the same as that for the AdS vacua (eqn. I167P provided we 
replace p in (fT67j) by p = A x QjA 2 P in (fTHTD and P - Q in (fWD by Q - P in (fTHTD . 
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For the de Sitter vacua, we use a less restrictive upper bound p = -^p < 100. For a manifold with 
N = 100 moduli we obtain: 

1) Ai = 3% , (188) 

2) A 2 = 31% . (189) 

In Figure [12] we present plots of log 10 (?7i3/2) for both cases as a function of P in the range where 
Peru < p < 100 for different values of Q — P. Recall that the smallest possible value of Q — P for de 

Log w (m V7 ) Log L0 (m V2 ) 




FIG. 12: log 10 (7713/2) as a function of P for case 1) - top plots and case 2) - bottom plots. The light grey area represents 
possible values of log 10 (m 3 /2) in the range where p cr %t < p < 100 consistent with the SUGRA approximation. The dark 
area indicates the region of interest where —16 < log 10 (m 3 / 2 ) < —14 such that 240 GeV < m 3 / 2 < 24TeV. The dark solid 
curve corresponds to log 10 (m3/2) when p = p C rit- The lower boundary of the light grey area represents the log 10 (m3/ 2 ) 
curve when p = 100. The dashed curve corresponds to p = 1.01. Top left: Case 1) when Q — P — 3. Top right: Case 1) 
when Q — P = 6. Bottom left: Case 2) when Q — P — 3. Bottom right: Case 2) when Q — P = 4. 

Sitter vacua is {C^ — P\min — 3. In this case the region where 240 GeV ^ ^3/2 — 24TeV, exists for 
3 < Q — P < 6 for the anisotropic case in fll88j> . and for 3 < Q — P < 4 for the "democratic" case in 
(HMD - 

D. Small Cosmological Constant implies Low-scale Supersymmetry in dS Vacua 

In this subsection we will study the distribution of SUSY breaking scales in the de Sitter vacua which 
as we showed earlier, can arise when the hidden sector has chiral matter. In particular we will see that 
the requirement of a small cosmological constant leads to a scale of SUSY breaking of 0(1 — 100) TeV. 
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In section VI, we saw that the minimum obtained is de Sitter if the discriminant of the quadratic 
polynomial with respect to <pQ in eqn. (|14ip is negative, while it is anti-de Sitter if the discriminant in 
(|14ip is positive. For 7713/2 ~ 0(1 — lOTeV) the magnitude of the vacuum energy in both cases can be 
estimated to be 

I Vol ~ m 2 p ml /2 ~ (10 10 GeV) 4 - (10 n GeV) 4 . (190) 

On the other hand, if the discriminant in (I141D vanishes, one obtains a vanishing cosmological constant 
(to leading order in the approximation). At present it is not known if there is a physical principle which 
imposes this condition. However, one can still use it as an observational constraint since the observed 
value of the cosmological constant is known to be extremely small. For instance, it could happen that 
the space of G2 manifolds scans the constants (Ai, P, Q, N) finely enough such that there exist vacua for 
which the vacuum energy is acceptably low. In particular, the constants A4, i = 1,2 which are determined 



by the threshold corrections have been shown to depend on integers 5 19]. In this work, we will assume 
that to be the case. A detailed computation to show this in a convincing manner is currently being 
attempted and will be reported in the future. It should also be kept in mind that a different mechanism 
for solving the cosmological constant, completely decoupled with particle physics, could exist. Such a 
mechanism, if present, would not affect any predictions for low energy particle physics. 
By setting the left hand side in (|144p to zero, we can then express 

^ [n \A 2 Pj 3(Q-P)-8' Uyij 

Of course, since the above constraint was obtained in the leading order, the vacuum energy is only 
zero in the leading order in our analytic expansion. The subleading contributions we neglected in (225) 
although smaller than the leading contributions, are still much larger than the observed value of the 
cosmological constant. However, one can in principle take into account all the subleading corrections and 
tune the ratio A1Q/A2P inside the logarithm to set the vacuum energy to a very small value compatible 
with the observations. As will be seen later, since the expression in the R.H.S. of (|19ip turns out to be 
large, the subleading corrections which affect the value of the cosmological constant will have little effect 
on the phenomenological quantities calculated by imposing the constraint to leading order. 

We would now like to analyze in detail the phenomenological implications of the solutions obtained 
by imposing (|19ip as a constraint. The most important phenomenological quantity in this regard is 
the gravitino mass as it sets the scale of all soft supersymmetry breaking parameters. We focus on the 
gravitino mass in this section. The soft supersymmetry breaking parameters will be discussed in section 

Em 

1. Gravitino mass with a small positive cosmological constant 
By substituting the constraint (I19ip into the gravitino mass formula (|187j) . we obtain: 
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where the meson vev is now given by: 



5 This is because the threshold corrections can be related to certain topological invariants of the associative three-cycle. 
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In this case, the moduli vevs are given by 
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From (|192p . one notes that the gravitino mass (when the cosmological constant is made tiny) is 
completely determined by the dual coxeter numbers of the hidden sector gauge groups A^andQ, the 
rational numbers a% (see ([2])) characterizing the volume of the G2 manifold and the integers iVj. 6 

The rationals a, are subject to the constraint 5_/i=i a i = It is reasonable to consider a "demo- 
choice for a,- 



cratic" choice for ctj, = 7/(3 A/") for all i = 1,N and also to take for simplicity all the integers N = 1. 
The integers Ni will generically be of 0(1); even if some of the Ni are unnaturally large, their individual 
contributions will be typically washed out as they are raised to powers that are much less than unity 
(see (|187p and the expression for aj for the democratic choice). In this case, after setting A2 = QC2, the 
gravitino mass formula is given by 
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and the moduli vevs are 
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From (|195p . the gravitino mass depends on just four constants - C2, P, Q and iV (the total number of 
moduli), determined by the topology of the manifold. It should be kept in mind that for the solution to 
exist, it is necessary that Q — P > 2 (see (|135|) ) . For the smallest possible value (Q — P) m in = 3, the 
expression for m 3 / 2 simplifies even further 
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Note that the dependence on A^ and P in (I197P is due solely to the volume Vx dependence on those 
parameters. The expression for the gravitino mass has a more transparent form if we don't substitute 
the expression for the volume (1186[) into A185 j> . For Q — P = 3 we obtain 
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P 8^ 3/2 ' 



-28 



Co ~ 514 TeV 



Co 



V 



3/2 



(198) 
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where the detailed dependence on aj, Ni, P and the number of moduli N is completely encoded inside 
the seven-dimensional volume Vx which appears to be the more relevant physical quantity. Furthermore, 
in the supergravity approximation the volume Vx > 1, which translates into an upper bound on the 
gravitino mass when Q — P = 3 



m 



3/2 



< 0(100 TeV) 



(199) 



6 From field theory computations [25( |. Ai = N c — Nf — P and A2 — Q (in a particular RG scheme), up to threshold 
corrections. We can therefore express Ai , A2 as A\ = P Ci and A2 — Q C2, where coefficients C\ and C2 depend only on the 
threshold corrections and are constant with respect to the moduli [igj | . In this case, the quantity ln(AiQ /A2P) = ln(Ci/C2) 
is fixed by imposing (|I9ip . 
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2. The Gravitino mass Distribution and its consequences 



The gravitino mass (|195j) depends on three integers: the two gauge group dual coxeter numbers N c 
, Q and the number of moduli N. 7 This will give us an idea about the distribution of the gravitino 
mass (which sets the superpartner masses) obtained after imposing the constraint (|19ip that the vacuum 
energy is acceptable. Asides from only considering the vacua within the supergravity approximation (ie 
Si > 1) we expect an upper bound on the dual coxeter numbers of the hidden sector gauge groups P and 
Q. Based on duality with the heterotic string, it seems reasonable to assume that they can be at least 
as large as 30 - the dual coxeter number of E%. Of course, values of P, Q larger than 30 cannot be ruled 
out, and here we assume an upper bound P < 100. Notice that from the constraint in (|19ip . the ratio 
(AiQ) / (A2P) can get exponentially large when Q — P = 3 and the values of P are small. Here we are 
going to completely relax the requirement on the upper bound of the ratio because as we will see, for 
generic manifolds with a large numer of moduli, P has to be large for the supergravity approximation 
to hold. The distribution can be constructed as follows. The three integers: P, Q — P and N are varied 
subject to (1135j) and the supergravity constraint Si > 1. For each point in the resulting two dimensional 
subspace, log 10 (7713/2) can De computed and rounded off to the closest integer value. One can then 
count how many times each integer value is encountered in the entire scan and plot the corresponding 
distribution. 
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FIG. 13: The gravitino mass distribution with the x-axis denoting the logarithm of the gravitino mass (to base 10). Left: 
Distribution corresponding to scan one in (|200p . Middle: Distribution corresponding to scan two in (|2U1 [I for which manifolds 
with the number of moduli TV < 50 were excluded from the scan. Right: Distribution corresponding to scan three in (|202[) 
for which manifolds with the number of moduli N < 100 were excluded from the scan. 



In the first three scans scan we cover a broad range of values by choosing P ma x = 100. Taking into 
account the SUGRA constraint (s, > 1), we have the following ranges of integers for the first scan: 



3 < P < 100 ; 3 < (Q - P) < 100 - P ; 2 < N < 



14 (P + (Q-P)) 
vr(3(Q-P)-8) 



(200) 



In the second scan we have excluded the small N region and considered only the manifolds with N > 50. 
Thus we have the following ranges of constants for the second scan: 



3<P<100; 3 < (Q-P) < 100 -P; 50 < A < 



14 (P + (Q-P)) 
tt(3(Q-P)-8) 



(201) 



In the third scan we have only considered manifolds with A > 100. Thus we have the following ranges 
of integers for the second scan: 



3 < P < 100 ; 3 < (Q - P) < 100 - P ; 100 < N < 



U(P + (Q-P)) 
vr (3(Q - P) - 8) 



(202) 



We can set C2 = 1 for the order of magnitude estimates we are doing here. 
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The first two distributions in Figure [13] clearly have several prominent peaks. Amazingly, in all three 
plots one of the peaks landed right in the m%/ 2 ~ 0(1 — 100) TeV range! The high scale peaks on the 
left plot appear to be around m 3 / 2 ~ 10 14 GeV and the GUT scales. However, for the middle plot the 
GUT scale peak almost disappears. Recall that the middle plot corresponds to scan two in (I20ip where 
we excluded all the manifolds for which the number of moduli N is less than 50. Therefore, the high 
scale peaks are largely dominated by contributions from the G 2 manifolds with a small number of moduli 
N < 50. As seen from the right plot, when G 2 manifolds with N < 100 are excluded from the scan, the 
peak at the m 3 / 2 ~ 10 14 GeV scale has all but disappeared, whereas the peak at m 3 / 2 ~ 0(1 — 100) TeV 
remains virtually unchanged. 

In Figure[l4]we included three more scans for which the upper bound on P was reduced to P ma x = 30. 
The fourth scan has the following ranges: 

14 (p _|_ (n — p)\ 

3<P<30; 3<(Q-P)<30-P; 2 < N < - V _ » ■ (203) 

In the fifth scan we again excluded the small N region an considered only the manifolds with N > 50 
and considered P ma x = 30: 

1 4 (p _i_ (n — p)\ 

3<P<30; 3< (Q-P) < 30 - P ; 50 < N < - p) _ 8) • (204) 

In the sixth scan we considered only the manifolds with N > 100 and P ma x = 30: 

14 (p _i_ (Q — p\\ 

3<P<30; 3< (Q-P) <30-P; 100 < N < \ -f g ■ (205) 

7r (3(Q - P) - 8) 

Again, in Figure [TJ] we notice that the 0(1 — 100) TeV peak narrows around m 3 / 2 ~ 0(100) TeV, as 




-3. -10 -15 -20 -25 -5 -10 -15 -20 -25 -5 -10 -15 -20 -25 



FIG. 14: The gravitino mass distribution with the x-axis denoting the logarithm of the gravitino mass (to base 10). Left: 
Distribution corresponding to scan four in (|203[) . Middle: Distribution corresponding to scan five in (|204[) for which manifolds 
with the number of moduli iV < 50 were excluded from the scan. Right: Distribution corresponding to scan six in (|205|l for 
which manifolds with the number of moduli TV < 100 were excluded from the scan. 



we exclude manifolds with small number of moduli. As the same time, the peaks at the high scale 
completely disappear for G 2 manifolds with N > 50. Finally, in Figure [151 we chose the smallest possible 
value Q — P = 3 and scanned integers P and N in the following ranges: 

14 (P + 3) 

3 < P < 200 ; 50 < N < — , (206) 

7T 

14(P + 3) 
3 < P < 100 ; 50 < N < — i '- , 

7T 

14 (P + 3) 
3<P<30; 50<iV< — 1 ; . 

7T 
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FIG. 15: The gravitino mass distribution with the x-axis denoting the logarithm of the gravitino mass (to base 10). Scans 
for the smallest possible choice (Q — P)min = 3. Left: Distribution corresponding to the scan with Pmax = 200. Middle: 
Distribution corresponding to the scan with Pmax = 100. Right: Distribution corresponding to the scan with Pmax = 30. 

In all three plots in Figure [T5l we see the same peak at m 3 / 2 ~ 0(1 — 100) TeV, which narrows around 
m 3/2 ~ 0(100) TeV as P ma x is decreased. 

Therefore, from the above distributions we conclude that the peak corresponding to ~ (1 — 

100) TeV is entirely due to the smallest possible value (Q — P) m i n = 3. This can be explained if we 
examine the gravitino mass formula in (|197H . In particular the constant factor e -28 ~ 10 -12 is most 
crucial in lowering the gravitino mass to the TeV scale. It is easy to trace the origin of this factor to the 
constraint (j!91|) . imposed by the requirement to have a zero cosmological constant (to leading order). 
When (|19ip is used along with the requirement Q — P = 3 we simply get 

P b (|f)=84. (207) 

When this is substituted into the gravitino mass (|187p . the corresponding suppression factor turns into 
the constant 



AtQ \ Q-p 
AoP 



e~ 28 . (208) 



Physically, this suppression factor corresponds to the hidden sector gaugino condensation scale (cubed). 
Recall that for an SU(Q) hidden sector gauge group, the scale of gaugino condensation is given by 



A g = m p e =m p e s5 Im/ . (209) 

The moduli vevs in (I194|) completely determine the gauge kinetic function. Taking Q — P = 3 we obtain 

N . n 

Im/ = y^N lSl = — ^. (210) 

i ✓ yj- 

i=l 

Substituting (|210p into (|209|) we obtain the following scale of gaugino condensation 

A g = m p e- 2S/3 « 2.15 x 10 14 GeV. (211) 

It is important to note that the expression in the R.H.S. of (|19ip is quite large (= 84, when Q — P = 3) in 
the leading order, and the quantity (A1Q/A2P) which is fixed by imposing the vacuum energy constraint 
is inside a logarithm. Therefore, even when one incorporates all the higher order corrections and tunes 
the ratio A1Q/A2P inside the logarithm to set the cosmological constant equal to the observed value, 
the constant on the R.H.S. (=84), crucial in obtaining the 0(100) TeV scale peak, is hardly affected. 
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FIG. 16: Plot of log 10 (m3/2) as a function of P for Q — P 
curve corresponds to N — 500. 



3. Short-dashed curve corresponds to N — 50. Long-dashed 



The dominance of the 0(1 — 100) TeV range also becomes clear from Figure [THl where log 10 (wi3/2) as 
a function of P for Q — P = 3 is plotted for a manifold with N = 50 moduli - short-dashed curve, and 
a manifold with N = 500 moduli - long-dashed curve. Indeed, even when we do not impose the SUGRA 
constraint, from the above plot we can see that the 0(1 — 100) TeV range is covered by a large swath 
of the vacuum space and it is not so surprising that the corresponding distribution peaks at that scale. 
This essentially follows from the formula for the gravitino mass in (|197p . 

An important point which should be emphasized is that for Q — P = 3, the gravitino mass dependence 
on P and N appears only through the volume Vx, as can be seen from (I198D . Thus, the distribution in 
Figure [15] directly correlates with the corresponding distribution of the stabilized volume of the seven- 
dimensional manifold Vx as a function of P and N. Therefore, it is the dominance of the vacua with a 
relatively small volume, which results in the peak at 0(100) TeV. 

Also note that in the above analysis we simply set the constant coefficient due to the threshold 
corrections 02 to unity. It would be interesting to get a handle on this quantity and include its variation 
into the gravitino mass distribution study. 

One could argue that even though Q — P = (Q — P)min = 3 gives a peak for the gravitino mass 
distribution at around 0(100) TeV scale, it seems plausible from a theoretical point of view to have 
many examples of gauge singularities in G2 manifolds such that Q — P > 3. However, by imposing the 
supergravity constraint that all moduli Sj are larger than unity (which is the regime in which the entire 
analysis is valid), one sees that having Q — P > 3 drastically reduces the upper bound on N compared 
to that for Q — P = 3 (see eqn ([194p ). Therefore, the peaks in the gravitino mass distribution obtained 
for —2 < log 10 (?7i3/2) < ~~ 5 in Figures [131 and [TH come from vacua with a small number of moduli as 
well as Q — P > 3, compatible with the analysis in the supergravity regime. Since it is presumably 
true that the number of G2 manifolds with the required gauge singularities which have a large number 
of moduli is much larger than those with a small number of moduli, it seems reasonable to expect that 
the peak of the gravitino mass distribution obtained at around 0(100) TeV scale is quite robust and is 
representative of the most generic class of G2 manifolds with the appropriate gauge singularities. Notice 
also that in the case of manifolds with a large number of moduli N > 100, because of the constraint 
from the supergravity approximation, the actual minimal value of P in the scans (|202p and (|205j) is quite 
large, i.e. P m i n = 20. Hence, if P > 20, from (|19ip the value of the ratio when Q — P = 3 is bounded 
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from above 

£ 67.7 (212) 

A 2 P 

On the other hand, for less generic manifolds with a small number of moduli the supergravity constraint 
allows P to be small. In this case in order to satisfy (|19ip . (AiQ)/(A2P) has to be exponentially large. 
Whether this is possible remains to be seen. However, because for the more generic case when N > 100 
the upper bound in (I212p is quite reasonable, our analysis remains robust. 

One could also contrast the results obtained above with those obtained for the Type IIB flux vacua. 
In Type IIB flux vacua, one has to independently tune both the gravitino mass to a TeV scale (if one 
requires low scale supersymmetry to solve the hierarchy problem) as well as the cosmological constant to 
its observed value. This is quite different to what we are finding here. Finally, we should emphasize that 
imposing the supergravity approximation was crucial in obtaining low scale SUSY breaking. Plausibly, 
the vacua which exist in the M theoretic, small volume regime will have a much higher SUSY breaking 
scale. However, such vacua presumably also have the incorrect electroweak scale (either zero or Mu). 



E. Including more than one flavor of quarks in the hidden sector 



In the previous analysis we assumed a single flavor for the quarks in the hidden sector, i.e Nf = 1. 
In order to check that the way we obtain a dS metastable minimum is robust and not dependent on a 
particular choice of chiral hidden sector matter spectrum, we would like to extend our analysis to include 
more than one flavor (but still with Nf < N c ) so that the meson fields are given by: 



2 Q a Qa 



1/2 



(213) 



where a, a = l,Nf. In the absence of a perturbative superpotential (which is guaranteed in the absence 
of fluxes), one has along the D-flat direction: 



Q = Q 



V2 



(<i>\ 



\ 



\ 



h N f 



Thus, the determinant appearing in (IMjl becomes 



det((f>%) = JJ (j> a , where <f> (T = <j)%. 



cr=l 



The nonperturbative superpotential and the Kahler potential are then given by 

Nf 



W = A X Y[<t> a a e iblf + A 2 e ib2f 
<j=i 

Nf 



(214) 



(215) 



(216) 



CT=1 
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where we again denoted b\ = 2tv/P, b 2 = 2n/Q, P = N c — Nf and a = —2/ P. After minimizing with 
respect to the axions, the scalar potential is given by 



V = 6 4 ^ y3 " [(61 A 1 n(^o)^- 6l " S - b 2 A 2 e- b ^«) 2 £ ciMf + 3(Ai W^e^ d - A 2 e~ b ^ d ) 2 

X (7=1 1=1 (7=1 

N f N f 

+3(z7- a)(A x Yl^ote-^ 3 - A 2 e- b ^ 3 ){b l A 1 \{(^) a a e- b ^ - b 2 A 2 e~ b ^ d ) (217) 

(7=1 (7=1 



7=1 " ' ' (7=1 

Instead of presenting a full analysis of this more general case we would simply like to check that we have 
a metastable dS vacuum, and that the main feature of the dS vacuum, namely the emergence of the TeV 
scale when the tree-level cosmological constant is set to zero, survives when Nf > 1. For this purpose we 
need to compute the scalar potential at the minimum with respect to the moduli s, as a function of the 
meson fields <$> a . 

The generalization of the equations minimizing the scalar potential is fairly straightforward. In par- 
ticular, in the limit when the size of the associative cycle Im/ = V- a is large, for A = 1 the generalization 
of the second equation in (|126p . which determines L\ + takes on the following form 

where we again defined (5 = ^IlSiC^o)? 6 "^" 62 ^" 3 ! x = (3 - 1, y = h(3 - b 2 and z = b\fi-b\. Thus, in 
the large three-cycle limit we again have /3 ~ &2/&1 = P/Q so that y — > and the leading order solution 
for -£/i,+ is again given by 

(218) 

In this case the moduli are stabilized at the same values given by (I132p . Since both the superpotential 
and the Kahler potential are completely symmetric with respect to the meson fields, it seems reasonable 
to expect that there is a vacuum where all <$> a are stabilized at the same value, i.e. ((po) a = 4>o f° r an 
a = 1, Nf. Using the solution for the moduli vevs (|132p and the above assumption we obtain the following 
expression for the potential at the extremum with respect to the moduli Sj as a function of 4>o 



V = N f ^ 



■ 2 ■ (2,1 ] 3 \ -., (u i s ~ 



•o + ( — - jf f ) <Po + 



1 " ' . (219) 



6g V 



By setting the discriminant of the biquadratic polynomial in the square brackets to zero we again 
obtain the leading order condition on the tree-level cosmological constant to vanish: 

3 8 28 , . 

Since the solutions for the moduli in the dS case correspond to branch b) where Q > P and A\Q > A 2 P, 
zero vacuum energy condition (|220p can be satisfied only when 
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Therefore, a vanishing tree-level cosmological constant in the leading order results in the following set of 
conditions: 



Pin 



AiQ\ _ 28(Q - P)N f 
A 2 P )~ 3(Q - P) - 8N f 



and (Q-P) 



(222) 



Recall that the key to obtaining the TeV scale gravitino mass was the exponential suppression factor 
e~ 28 when Q — P = (Q — P)min = 3, related to the scale of gaugino condensation. In the present case, 
up a factor of order one, we have 



in 



3/2 



m p ( AiQ\ 
Vl 12 \MP) 
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3/2 



3(Q-P)-8Nf 



'X x ' v x 
Consider a few examples where Nf > 1. From (|222j) we have the following set 
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N f 



N f 
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4, (Q-P)r 
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Nf = 4, Q-P = 12, 
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(224) 
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Remarkably, in all but one cases listed above we obtain the same suppression factor e -28 ~ 7 x 10~ 13 
which was the reason for the peak at m 3 / 2 ~ C(l — 100)TeV! Note that the only example which did not 
fall into this range was the third case for which the condition on the cosmological constant to vanish was 
Pin (^iS) = 1232, which is too unrealistic anyway, as it requires either extremely large dual coxeter 

numbers for the gauge groups N c , Q ~ 0(1000) or an exponentially large ratio inside the logarithm. On 
a similar note, as can be seen from the third entry in each line in (|224p . increasing the number of flavors 

Nf even further would again require either P,Q> 300 or an extremely large ratio (^^73^ j which appears 
inside the logarithm. Therefore, limiting our analysis to the cases with Nf < 5 seems quite reasonable. 

Recall that for Nf = 1 the TeV scale appeared for the minimum value (Q — P) m in = 3 whereas the 
vacua corresponding to the higher values of Q — P generally failed to satisfy the SUGRA constraint for 
more generic G 2 manifolds with a large number of moduli. For this reason, considering larger values of 
Q — P for the examples listed above is probably unnecessary. Hence, for more than one flavor of quarks, 
we only need to take 



Q-P = W 



f 



(225) 



Thus, given that the assumptions we made in the beginning of this subsection are reasonable, it appears 
that the connection of the TeV scale SUSY breaking to the requirement that the tree-level vacuum energy 
is very small is a fairly robust feature of these vacua, independent of the number of flavors. 



F. Including matter in both hidden sectors 



In the previous analysis we tried to be minimalistic and included chiral matter in only one of the 
hidden sectors. Due to this asymmetry, we obtained two types of solutions - a supersymmetric AdS 
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extremum when P > Q corresponding to branch a) and a dS minimum when Q > P (when condition 
(|144l holds), corresponding to branch b). Using this result it is then fairly straightforward to figure out 
what happens when both hidden sectors produce F-terms due to chiral matter. For the sake of simplicity, 
we will again consider the case when Nf = 1 in both hidden sectors. In this case, the Kahler potential is 
given by 

K = -3 ln(4vr 1 / 3 V x ) + # + ip$ . (226) 

After minimizing with respect to the axions, the non-perturbative superpotential (up to a phase) is given 
by 

W = -A l( j) ai e~T- Im/ + A 2 tfi a2 e~^ lmf , (227) 

where a\ = —2/P and a 2 = —2/Q. We will now check to see if it is still possible to obtain SUSY extrema 
when both hidden sectors have chiral matter. Setting the moduli .F-terms to zero we obtain 

3PQJ3-1 

Vk = v = ~^rw^^ (228) 

where (3 = ^~ "o"-* 1 ™^. At the same time, setting the matter F-terms to zero results in the 

following conditions: 

§ + l)/3-l = 0. (229) 

- 71 + P ~ 1 = • ( 23 °) 

Expressing (3 from (I228p and substituting it into (|229H230j) and using the definitions for a\ and a 2 , we 
obtain the following expressions for the meson field vevs: 

^ = 2 ± f^_ (231, 
rf-£±*2£M (232) 

Since v as well as both </>q and ipQ are positive definite, we have the following two possibilities: 

a) P >Q : = and ^ 0, (233) 

b) P <Q: F^ / and F$ = 0. 

Thus, when both hidden sectors have chiral matter, supersymmetric extrema are absent. Instead, when 
condition (|144p holds (for branch a) we simply swap P and Q, A± and A 2 in (I144p ). for each branch we 
obtain a dS vacuum where only one of the matter F-terms is non-zero. Keep in mind that although in 
the above analysis we used condition (|228h obtained by setting the moduli F-terms to zero, even in the 
dS case when the moduli F-terms are non-zero, one of the two mesons will be stabilized at a value such 
that the corresponding matter F-term is zero. The zero F-term has no effect on the analysis of the dS 



53 



solution and apart from replacing a with j3 denned above, the same solution obtained previously for the 
dS vacuum applies. In this case, the only difference will be in the meson field vevs: 



(234) 



Q-p oinf^QV™ Q-p V Q-p P UMg)\2 V Q-p 



A 2 P J * LL1 \A 2 P 

Therefore, the dS solution obtained for the minimal case when only one of the hidden sectors has chiral 
matter does not change even when we include chiral matter in both hidden sectors. 

G. Height of the Potential barrier 

For simplicity, we first compute the height of the potential barrier for the case with a pure SYM 
hidden sector. The two solutions for A = 1 in eq.(50) and eq.(51) exist for any number of moduli and 
therefore the analysis below extends to the general case with an arbitrary number of moduli. The solution 
in eq.(51) corresponds to the dS maximum which determines the height of the barrier. Using eq.(44), 
from eq.(51) we have 

4 = **(*±£E£\. (235, 



7 y \ w 

Therefore, we can express the volume of the associative cycle Vol{Q) = V ■ a as 



x 



3 + V9 - 7w 



u-a= , (236) 

y V w J 

where we used the fact that v\ = v in (|235j) is independent of k and X^=i a * = |- Using the definitions 
of x, y, z, and w in eq.(37) in terms of a, from (|236j) we can solve for a 

P 2 (lQ 2 + 12ttQz7 • a + 4tt 2 (V ■ a) 2 ) 

a = ) -( . (237) 

q2 l 7P 2 + ■ a + 4vr 2 (V ■ a) 2 J 

From eq.(36) we can express v ■ a as 

^ ^ 1 PQ (A 2 \ 

v ■ a = — — — In —a . (238) 

In the limit when the volume of the associative cycle V ■ a is large we can solve (1237)) and (|238p order by 
order to obtain 

■■Wi, 6(P-Q) 2 {29P-1Q)(P-Qf 



together with 



, p g' n (w) 3(P-Q) ll(P-Q)\P + Q) 
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The moduli vevs at the barrier are then given by 




3(P-Q) , 11(P- Q) 2 (P + Q) 



MP-Q) ".ln(^) + 2nPQln 2 (^ 



(241) 



In the leading order 
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7Ni 



2it(P - Q) 



(242) 



and 
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P2 

Q 2 ' 



(243) 



Using (1242p and (|243p . the value of the potential at the barrier in the leading order is given by 
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(244) 



Recall that the value of the gravitino mass at the SUSY AdS extremum in the leading order is given by 
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(245) 



Therefore, we can express the value of the potential at the barrier in the leading order as 



2 2 
m p m 3/2 



7{P 2 - Q 2 ) 2 + PQ In 



7Q 2 (P-Q) 2 \Q 
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2P 
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7 (P 2 - Q 2 ) + PQ In 



A 2 P< 
AiQ 2 



(246) 



Note that the above expression is independent of the number of moduli and parameters and iVj. This 
formula is very accurate when compared to the exact numerical values. For example, for the values of 
the parameters in eq.(54), the numerically obtained result is 



Vb ~ 51.55 x m 2 m 2 / 2 , 
whereas from the leading order expression in (|246p we get 

V b w 49.92 x m 2 m^ /2 . 



(247) 
(248) 



In Figure 1171 below we plotted the value of the scalar potential at the barrier in units of 

Tri^TTi^^ as a 

function of P for P — Q = 3 and P — Q = 1 and two different values of parameter p = A^ . 

For the case of vacua with charged matter in the hidden sector, the situation is more complicated as 
the potential depends on an additional field (meson) . The height of the potential barrier changes as one 
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FIG. 17: Plots of the potential at the barrier in units of m^m^ as a function of P. Solid line corresponds to -j^J = 100 
while the dashed line corresponds to = 10. Left: Plot for P = Q + 3. Right: Plot for P = Q + 1. 




FIG. 18: Potential in units of m\i 2 m^, plotted as a function of the meson field and the moduli along the slice si = S2, for 
a manifold with two moduli and the microscopic constants as in (|136[) . 

moves along the meson and the moduli directions. To illustrate this with an example, we have done a 
numerical analysis for a manifold with two moduli and one meson field with the microscopic constants 
as in (|136p . Figure [18] shows a three-dimensional plot of the potential as a function of the moduli along 
the slice s\ = S2, and the meson field <j). In figure [T9l in the left plot the potential is plotted along the 
slice s\ = S2 with the meson field equal to its value at the minimum (|137p . while in the right plot the 
potential is plotted along the slice s\ = S2 with the meson field equal to a value such that the height of 
the potential barrier is at its minimum. 

VIII. PHENOMENOLOGY 

In this section, we will begin the analysis of more detailed particle physics features of the vacua, 
with emphasis on the soft supersymmetry breaking parameters, since we are particularly interested in 
predicting collider physics observables that will be measured at the LHC. 

The low energy physics observables are determined by the kahler potential, superpotential and the 
gauge kinetic function of the effective N = 1, d = 4 supergravity. The gauge kinetic function (/) has 
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FIG. 19: Left: Potential in units of m 3 / 2 m p along the slice si = S2 for a manifold with two moduli with the meson field 
equal to its value at the minimum of the potential (as in fl37j)). Right: Potential in units of m^ 2 m p along the slice si = S2 
for a manifold with two moduli with the meson field <f> = 0.102 (this is such that the height of the potential barrier is at its 
minimum). The microscopic constants for both cases are the same as in (|136[l . 



already been discussed. The Kahler potential and the superpotential can be written in general as follows: 



K = K{s h< j> h ,4> h ) + K &f3 (si)$^ + Z aP ( Si ,cf> h )$ a $P + ... 
W = W(zi) + //$ a $£ + Yt tfh $ a $ f} & + ... 



(249) 



where <3? a are the visible sector chiral mater fields, K^p is their Kahler metric and Y' a ^ are their unnor- 
malized Yukawa couplings. <p h denote the hidden sector matter fields. The first terms in K and W depend 
only on the bulk moduli and have been already studied earlier. In general there can be a mass term (//) 
in the superpotential, but as explained in [2(i], natural discrete symmetries can exist which forbid it, in 
order to solve the doublet-triplet splitting problem. The quantity Z a p in the Kahler potential will be 
important for generating an effective fi term as we will see later. 

Since the vacua have low scale super symmetry, the effective lagrangian must be equivalent to the 
MSSM plus couplings involving possibly additional fields beyond the MSSM. For simplicity in this section 
we will assume an observable sector which is precisely the MSSM, although it would also be interesting 
to consider natural M theoretic extensions. The MSSM lagrangian is characterized by the Yukawa and 
gauge couplings of the standard model and the soft supersymmetry breaking couplings. These are the 
scalar squared masses m?, the trilinear couplings Aijk, the fi and B/i mass parameters and the gaugino 
masses. In M theory all of these couplings become functions of the various constants (Aj, N, P,Q, iVfe) 
which are determined by the particular GVmanifold X. In addition, because we are now discussing the 
observable sector, we have to explain the origin of observable sector gauge, Yukawa and other couplings 
in M theory. As we have already explained, all gauge couplings are integer linear combinations of 
the N moduli, the N integers determining the homology class of the three dimensional subspace of X 
which supports that particular gauge group. Furthermore, the entire superpotential is generated by 
membrane instantons, as we have already discussed. Therefore mass terms and Yukawa couplings in 
the superpotential are also determined by integer linear combinations of the moduli fields. Hence, in 
addition to the constants (A%, N, P, Q, N^) which determine the moduli potential, additional integers 
enter in determining the observable sector superpotential. Generically, though, we do not expect these 
integers to be large in the basis that the moduli Kahler metric is given by (6). 

We determine the values of the soft SUSY breaking couplings at M un if in the standard way : The 
moduli fields, hidden sector matter fields as well as their auxiliary fields are replaced by their vevs in the 



M = l,d = 4 SUGRA lagrangian. One then takes the flat limit M, 



a global SUSY lagrangian with soft SUSY breaking terms 23] . Unfortunately, in M theory the matter 



oo with m 3 /2 fixed. 



This gives 
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Kahler potential is difficult to compute. This leads to theoretical uncertainties in the calculation of the 
scalar masses and A, B and \x parameters. Fortunately though we are able to calculate the gaugino 
masses. Our main phenomenological result is that the tree level gaugino masses are suppressed relative 
to the gravitino mass. After explaining this, we will go on to discuss the other soft terms in a certain, 
calculable limit. 



A. Suppression of Gaugino masses 



Grand Unification is particularly natural in G 2 vacua of M theory [19|] . This implies that the gaugino 
masses at tree level (at the unification scale) are universal, i.e. the gauginos of the three SM gauge groups 
have the same mass. In order to compute the SM sector gaugino mass scale at tree- level we need the 
Standard Model gauge kinetic function, f sm . In general this will be an integer linear combination of the 
moduli, with integers Nf m , which is linearly independent of the hidden sector gauge kinetic function in 
general. The expression for the tree-level MSSM gaugino masses in general N = 1, d = 4 SUGRA is 
given by: 



K/2 f>nmp_ a f 
c - £i 1 m u n J s 

2i Im/ S . 



1/2 = m P „. T _ t , (250) 



Note that the gauge kinetic function is independent of the hidden sector matter fields. Therefore, the 
large hidden sector matter F term responsible for the dS minimum does not contribute to the gaugino 
masses at tree level. We will now proceed to evaluate this expression explicitly both for the AdS and dS 
vacua. We will find that generically, the gaugino masses are suppressed relative to the gravitino mass. 



1. Gaugino masses in AdS Vacua 

Choosing the hidden sector to be pure SYM with gauge groups SU(P) and SU(Q), the normalized 
gaugino mass in these compactifications can be expressed as 



M 



m„ e 



1/2 



80F^J /2 



P 



Q 



2-7T T £ 27T 

+ A ie -^ lraf - A 2 e~Q 



Im/ 



(251) 



where is the phase of the superpotential W. In the leading order, the last two terms in the brackets 
can be combined as 



A ie -^ lm f -A 2 e 



A 7 



P-Q 
Q 



A 2 P 
AiQ 



P-Q 



(252) 



On the other hand, the two terms in the round brackets coming from the partial derivative of the 
superpotential cancel in the leading order. Therefore we need to take into account the first subleading 
order contribution (|92ft . In this order, we obtain 8 : 



P 



Ao -2^, 

— e Q 

Q 



±A 2 
2vr 



P-Q 
Q 



A 2 P 
A±Q 



p=Q B 



A 



(253) 



Recall that Im(f)^ c) = T A (c) (see dJJl). 



58 



From (|252p and (|253j) we notice that the absolute value of gaugino mass can now be conveniently expressed 
in terms of the gravitino mass (for a given value of A and c, as discussed in previous sections) as 



\ Ai- 



de) 

1/2 U 



2 e£i n 4% (4% + 3/2) (atv^ 



-(c) 



1 ^( c ) 



(254) 



where we also used ([79)) and (1851) . Finally, using ([7311 and the first equation in (j72l) . after some algebra 
we arrive at the following expression for the gaugino mass: 



\M- 



i(c) 
1/2 U 



3 



TX ; + 1 



q-A 



%f X (^3/2)a . 



(255) 



where we have introduced a new quantity 



such that the range of possible values for q is 



1 < q < 1 



(256) 



(257) 



Note that the general formula (1255[) which relates the gravitino and gaugino masses is completely 
independent of the number of moduli. 

When all rrik have the same sign the gaugino mass in (|255j) automatically vanishes. This is expected 
since the solution when A = ±1 is the SUSY extremum. In Figure [20] we have plotted absolute values of 



{Mx/2) a and (Mi/ 2 )a' as functions of q. 
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FIG. 20: Absolute values of (M 1/2 )^ l; -left and (M 1/2 )a ) -ri ght in units of gravitino mass as functions of q. As parameter 
A varies over < A < 1/7 - on the left and < A < 1 - on the right and the whole light grey region is covered. 
Left plot: A = 0- long dashed line, A = 1/9 - solid line, A — 5/36 - short dashed line. 
Right plot: A — 0- long dashed line, A — 0.5 - solid line, A = 0.95 - short dashed line. 



For a significant fraction of the space in both plots we have (Mi/ 2 )V < 0.2 (m 3/ / 2 )V , so the gaugino 
masses are typically suppressed compared to the gravitino mass for these AdS vacua. Note also that the 
suppression factor in (|255j) is independent of the gravitino mass. This result is different from the small 
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hierarchy between Mi/ 2 and m 3 / 2 in the Type IIB flux vacua ll|], where the gaugino mass is generically 
suppressed by In (7743/2). 

For the special case A = 0, system (|72p yields two solutions with positive moduli. Therefore, there 
will be two different values for the gaugino mass corresponding to these solutions. After some algebra we 
obtain 



\M , l (1 ' 2) 

\ 1V1 1/2\Q 



17 



9± 



-9+V17 



x ("13/2)0 



(1.2) 



0.22 



(^3/2)0 



(1,2) 



(258) 



Again, this relation is valid for any AdS vacuum with broken SUSY with A = with an arbitrary number 
of moduli. 

To check the accuracy of the approximate gaugino mass formula we again try the special case with 
two moduli ai = a 2 = 7/6 with the same choice of the constants as in (j67[) and the integer combination 
for the Standard Model gauge kinetic function {Nf m = 2, iV| m = 1}- In this case equation (I258|) for the 
absolute value of Af 1 / 2 yields: 



(M l/2 )V> = 164.4 GeV, (M 1/2 )^ = 95 GeV , 
whereas the exact values computed numerically for the same choice of constants are: 

(M 1/2 )o ] = 165.4 GeV, (M 1/2 ) 2) = 97 GeV . (260) 
This demonstrates a high degree of accuracy of our approximation, similar to that for the gravitino mass. 



*(2) 



(259) 



2. Gaugino masses in dS Vacua 

From the formula for the gaugino mass in (|250fl . the gaugino mass for the dS vacua in general can be 
expressed as 
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; L 1}+ + 1 m 3/2 , 



(261) 



where in the second equality we used (|122p and the fact that for these vacua vi = v for all % = 1,N, 
independent of i. Also, by including the minus sign we took into account that m 3 ^ 2 = e K / 2 \x\ A2e~ b2U ' a 
but x < 0, since Q — P > 3. From (|127p we can find L\ + including the first subleading contribution 



~ 3 3aabiy ( ad 

Li} 4. ~ 1 ( — „ h 1 

i,+ 2 14x5 



(262) 



For x, y and z in ()262p we use the definitions (|114p and substitute the approximate result f)133|) for a. 
Then after substituting (|262p into (I26ip and assuming that Q — P ~ in the limit when P is large 

the approximate tree level MSSM gaugino mass is given by 
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(263) 
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where we use ()134j) to substitute for 4>q. It is important to note some features of the above equation. First, 
equation (|263p is completely independent of the choice of integers Nf m for the Standard Model gauge 
kinetic function as well as the integers Ni for the hidden sector. Second, it is independent of the number of 
moduli N and moreover, it is also independent of the particular details of the internal manifold described 
by the rational numbers appearing in the Kahler potential ([2]). These properties imply that relation 
(|263l) is universal for all G2 holonomy compactifications consistent with our aproximations, independent 
of many internal details of the manifold. Furthermore, the denominator P \n.(AiQ / A 2 P) turns out to 
be always greater than unity for choices of microscopic parameters consistent with all constraints such 
as the supergravity regime constraint and the dS minimum constraint. In fact for reasonable choices 
of parameters, it is typically of O(10 — 100). Since the expression in the round brackets in (1263P is 
slowly varying and for the range under consideration is of 0(1), we see that gaugino masses are always 
suppressed relative to the gravitino for these dS vacua. 

After one imposes the constraint equation (|19ip ( (|207p when Q — P = 3) to make the cosmological 
constant very small, one can get rid of one of the constants in (|263p . and further simplify the expression 
for the universal tree level gaugino mass parameters for the dS vacuum with a very small cosmological 
constant: 

M 1/2 » 1 + —5 + — -= xm 3/2 = - e - irtw V xm 3 / 2 « -e~ l ™ 0.024 xm 3/2 . (264) 

1/2 84 V 3^ 84$ J 3 ' 2 34356 3/2 3/2 V ; 

As in the more general case (eqn. (|263p ). the tree level gaugino mass is suppressed compared to the 
gravitino mass, and the suppression factor can also be predicted. 

One would like to understand the physical origin of the suppression of the gaugino masses at tree level, 
especially for the dS vacua which are phenomenologically relevant. As mentioned earlier, since the matter 
F term does not contribute to the gaugino masses, the gaugino masses can only get contributions from 
the moduli F terms, which, as explained in the last paragraph in section fVI D|. vanish in the leading order 
of our approximation. The first subleading contribution is suppressed by the inverse power of the volume 
of the associative three-cycle, causing the gaugino masses to be suppressed relative to the gravitino. Since 
the inverse volume of this three-cycle is essentially cxhidden ~ 

the hidden sector gauge coupling in the UV - 
the suppression is due to the fact that the hidden sector is asymptotically free. In large volume type IIB 
compactifications, the moduli F terms also vanish in the leading order, leading to suppressed gaugino 



masses as well ll|]. However, in contrast to our case, there the subleading contribution is suppressed 
by the inverse power of the volume of the compactification manifold. Note that a large associative cycle 
on a G 2 manifold does not translate into a large volume compactification manifold. Thus, unlike large 
volume type IIB compactifications, these M theory vacua have a much higher compactification scale and 
hence are consistent with standard gauge coupling unification. 



B. Other parameters and Flavor issues 

The trilinears, scalars, anomaly mediated contributions to gaugino masses and the Bfj, parameter 
depend more on the microscopic details of the theory - the Yukawa couplings, the // parameter and the 
Kahler metric for visible sector matter fields. The flavor structure of the Yukawa matrices as well as that 
of the Kahler metric for matter fields is crucial for estimating flavor changing effects. We will comment 
on these at appropriate places. 

The (un-normalized) Yukawa couplings in these vacua arise from membrane instantons which connect 
singularities where chiral superfields are supported (if some singularities coincide, there could also be 
order one contributions). They are given by: 

^ = ^7^"^^ ( 265 ) 
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where C aj a 7 is an 0(1) constant and If are integers. 

The moduli dependence of the matter Kahler metric is notoriously difficult to compute in generic 
string and M theory vacua, and the vacua under study here are no exception. The best we can do here 
is to consider the Type IIA limit of these vacua. The matter Kahler metric has been computed in type 



IIA intersecting D6-brane vacua on toroidal orientifolds [28| building on earlier work 29]. Since chiral 



fermions livin g at intersections of D6-branes lift to chiral fermions supported at conical singularities in 
M theory 15|, 3(3] , we will simply uplift the IIA calculation to M theory. Thus the results of this section 
are strictly only valid in the Type IIA limit. 

Lifting the Type IIA result to M theory, one gets (see Appendix for details): 

K * - *<# 11 {-Twr) 

tan(7r#f) = cf(si) 1 ; cf = constant; I = rational number of 0(1). (266) 

In the type IIA toroidal orientifolds, the underlying symmetries always allow us to have a diagonal 
Kahler metric (28|. We have assumed for simplicity the Kahler metric to be diagonal in the analysis 
below. Now, we will write down the general expressions for the physical Yukawa couplings and the soft 
parameters - the trilinears and the scalars and then estimate these in M theory compactifications. The 
u and B/i parameters will be discussed in section IVIIIEI 

The Kahler potential for the chiral matter fields is non-canonical for any compactification in general. 
In determining physical implications however, it is much simpler to work in a basis with a canonical Kahler 
potential. So, to canonically normalize the matter field Kahler potential, we introduce the normalization 
matrix Q : 

$^Q-$, s.t. Q)KQ = 1. (267) 

The Qs are themselves only defined up to a unitary transformation, i.e. Q' = Q ■ U is also an allowed 
normalization matrix if IA is unitary. If the Kahler metric is already diagonal {K^p = K a 5a/3), the nor- 
malization matrix can be simplified : Q €xsi = (i^,) -1 / 2 ^. The normalized (physical) Yukawa couplings 



arc 



27j: 



Y al 3y = e X / 2 -^ w Q Q > Q Q W <2y 7 (268) 

It was shown in j^l] that in the class of theories with a hierarchical structure of the un-normalized 
Yukawa couplings (in the superpotential), the Kahler corrections to both masses and mixing angles 
of the SM particles are subdominant. Since in these compactifications, it is very natural to obtain a 
hierarchical structure of the un-normalized Yukawa couplings due to their exponential dependence on the 
various moduli and also because of some possible family symmetries, therefore one expects the effects 
of the Kahler corrections which are less under control, to be subdominant. The expressions for the 



un-normalized trilinears and scalar masses are given by 271 ]: 



m % = {ml /2 + V )K^-e k F^(d^d n k^-d m k^k^d n K- 5p )F n (269) 

Kfr = ^-z K/2 F m [K m Y' aM + d m Y^ - {K & -p d n K pa Y^ + a~ 1 + a~ 0)] 

The normalized scalar masses and trilinears are thus given by: 

ml/3 = (Q ] ■ m' 2 ■ Q)^ (270) 
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Let us discuss the implications for the soft terms, beginning with the anomaly mediated corrections 
to gaugino masses. 

Anomaly mediated contributions to Gaugino Masses 

We saw in section IVIH Al that the gauginos are generically suppressed relative to the gravitino. Since 
anomaly mediated gaugino masses are also suppressed relative to the gravitino (by a loop factor), they 
are non-negligible compared to the tree level contributions and have to be taken into account. Also, 
since anomaly mediated contributions for the three gauge groups are non-universal, these introduce 
non-universality in the gaugino masses at the unification scale. 

The general expression for the anomaly mediated contributions is given by (32| : 

(M)r = _A_ { _ {3Ca -Y,C2)e k ' 2 W* + (C - J2 C2)e k ' 2 F m K m + 2 ]T(C a a e^ 2 F m d m ln(£ a ))] 

a a. a. 

(271) 

where C a and C" are the casimir invariants of the a th gauge group and a runs over the number of fields 
charged under the a th gauge group. For a given spectrum such as that of the MSSM, C a and C™ are 
known. 

We first compute the F-term contributions 

e k ' 2 F i K i = ye-*w (i lt+ + ^jx m 3/2 (272) 



Then, equation (|271|) gives 



(M)r = _ e -i W ^[_ (3Ca _ ^ C:) + ^(C a - ^ C«) + ^] (273) 

+ (Ca-Y, C ") (f + 4) ~ \ + ~) E C " E i (* ^ Sin ( 2 ^))] X m 3/2 . 



where 



and we have defined the quantity: 
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{oigut)~ 1 = E s i N i m > (274) 



i=l 



ffl^, (275) 

where Of implicitly depends on the moduli. However, it is much simpler to keep the dependence as a 
function of 9f, as is explained in the Appendix. Depending on the values of the Casimir invariants C a 
and C % a for the three gauge groups, the anomaly mediated contribution can either add to or cancel the 
tree level contributions. Here we also took into account that TO3/2 = e K / 2 \x\ A2e~ b2U ' a but x < 0, since 
Q — P > 3. Using the expression for £14. in (|129p along with the definitions of x, y and z in (|114j) in 
terms of a in (|133j) and assuming that Q — P ~ 0(1), in the limit when P is large we obtain 

k - - -I + [ l + wh>* + hp*(0> 1 ■ (276> 
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Using (|143p and substituting for v from f)132|) into (|143j) together with a = —2/P we can express 

(277> 

Substituting (l276|) and ([277) into (l273j) we obtain 

(M)r _ _ e -,™^,_ (3a _ E< ~, + (, + _L_ + ) (278) 



x m 



3/2 • 



Note that these M theory vacua do not have a no-scale structure. Therefore, the anomaly mediated 
gaugino masses are only suppressed by loop effects, in contrast to the type IIB compactifications, which 
exhibit a no-scale structure in the leading order , leading to an additional suppression of the anomaly 
mediated gaugino masses. 

As before, when one imposes the constraint (|207p . the anomaly mediated gaugino mass contribution 
can be simplified further and is given by: 

(M)r = _aaur^ 1 ,_ (3Co _ E C) + ^055 + U374^ (Ca _ ^ ^ 



( M )r « - aG ^^ 7M/) [-(3Ca + 1-6556 (C a -^C) (279) 

- 0.048^ C a a £ i- (I C sin(2vr^))] x m 3/2 . 

a i 

From the left plot in Figure IA 21 of the Appendix we note that \^{l?P? sin(27T0f))| < 0.5. In a generic 
case, we expect that parameters Of are all different and, as a result, the terms appearing inside the 
corresponding sum over i partially cancel each other. Thus, in a typical case we expect that 



< 1 . (280) 



Neglecting the corresponding contribution in (|279p . taking ogut = 1/25, and substituting the Casimirs 
for an MSSM spectrum, we obtain the following values in the leading order, up to an overall phase e~' nw : 

(M)^i) « 0.01377 x mg/2, (M)^ (2) « 0.02317 x m 3/2 , (M)^ (3) « 0.02536 x m 3/2 . (281) 

Finally, combining the tree- level (|264p plus anomaly mediated (|28ip contributions, we obtain the following 
non-universal gaugino masses at the unification scale: 

Mi Kt -10.24 x 10~ 3 m 3/2 , M 2 rs -0.84 x 10~ 3 m 3/2 , M 3 +1.35 x 10~ 3 m 3/2 . (282) 

We immediately notice remarkable cancellations for M 2 and M 3 between the tree-level and the anomaly 
mediated contributions. Recall that since the distribution of ?n 3/ / 2 peaked at m 3 / 2 ~ 0(100) TeV, the 
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possible range of gaugino masses is in the desirable range m 1/ / 2 ~ 0(0.1 — 1) TeV. One of the consequences 
of these cancellations is a comparatively lighter gluino. Furthermore, since M 2 is a lot smaller than Mi, 
assuming R parity conservation, the neutralino LSP is expected to be wino-like. This is confirmed by 
explicitly RG evolving the gaugino masses to low scales, at least for the case when the cosmological 
constant is tuned to be very small. 

One should be extremely cautious however, since the predictive expressions above are only true if 
(|280p is satisfied. The extra contribution neglected in the above estimates is given by: 

A a = 0.15 x 10- 3 E C a V sin ( 2 ^)) x ™3/2 • (283) 

a i 

Because of the large cancellations between the tree- level and anomaly mediated contributions, it may 
happen that these corrections become important in a relatively small region of the overall parameter 
space, leading to a deviation from the above result thereby altering the pattern of gaugino masses. Further 
corrections may also come from varying ccgut as well as taking into account subleading corrections to 
the condition for the cosmological constant to be very small. In section IVHI C[ the effects of subleading 
corrections to the very small cosmological constant condition on the gaugino masses will be analyzed. A 
thorough study of these issues will be done in [131 ]. 



Trilinears 



The normalized trilinear can be written as : 



Aafr = ^ K ' 2 {KaKpK^)- 1 ' 2 (E e K ' 2 F m [K m + d m Y' a ^ - d m ln(K a k,K,)]\ (284) 



E eK/2pm [Km + d m ln(y^ 7 ) - 8 m HKaKpK^ 



As stated earlier, the subscripts {a, f3, 7} stand for the visible chiral matter fields. For example, a can be 
the left-handed up quark doublet, /3 can be the right-handed up quark singlet and 7 can be the up-type 
higgs doublet. Our present understanding of the microscopic details of these constructions does not allow 
us to compute the three individual trilinear parameters - corresponding to the up-type Yukawa, the 
down-type Yukawa and the lepton Yukawa matrices, explicitly. One can only estimate the rough overall 
scale of the trilinears. 

We see from ()284p that the normalized trilinears are proportional to the Yukawas since the Kahler 
metric is diagonal. If instead the off-diagonal entries in the kahler metric are small but non-zero, it would 
lead to a slight deviation from the proportionality of the trilinears to the Yukawa couplings. In most 
phenomenological analyses, the trilinears A are taken to be proportional to the Yukawas and the reduced 
trilinear couplings A a p y = A a01 /Y a01 are used. We expect this to be true in these compactifications 
from above. If the Yukawa couplings are those of the Standard Model, then from (|284p the normalized 
reduced trilinear coupling A a01 for de Sitter vacua in general is given by 



A a01 = e~ ilw I 1 + — 2 ^ |9 + r I Ul + 7 r [7 + 2 In 

1 (Q-P)<t>l d&Phi(4£\ Pin 4-^ 
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If we then use 
the trilinears 



together with (|266|) and Qa/3 = (K a ) 1//2 <5 a /3) we obtain the following expression for 
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Imposing the constraint equation (|207p on the expression above, the reduced trilinears for a dS vacuum 
with a tiny cosmological constant are simplified to: 



A 



a/3-y 



69 + 22^3 139 + 396^3 



-6 In 



7T 



72 

-E 



34356 



[7 + 21n 



a 



a/3"/ 



Y 



a/3"/ 



llnf^ 1 "^ 

2 V r(flf) 



1 

2^' 



Z V, ? sin(27r0f ) ^ + a -»• /3 + a -► 7 )]) x m 



l 3/2 



^4a/?7 



(1.4876 + 0.024 [10.45 + 2 In 



a 



a/3"/ 



Y 



a/?7 



7 In 



14(P + 3) 



N 



(287) 



2 V r(0?) 



sin(2vr^) f + «-/?.. . • j ,„« 2 . 



We see that compared to the gauginos, the trilinears depend on more constants. The quantity 
j 2 hi ( T< r(ef) ~ ) ~~ Sr(^f sin(27r0")| is of 0(1). Therefore, in a generic situation, we expect the terms 

inside the sum in |^hi ( ~T(gy^ ' ) ~~ ^(^f sin(27r#f)j to partially cancel each other and give an 
overall contribution much smaller than the first three terms inside the square brackets. Then, for known 
values of the physical Yukawa couplings and reasonable values of P and N, the trilinears generically turn 
out to slightly larger than m 3 / 2 . 

Scalar Masses 



For an (almost) diagonal Kahler metric, the normalized scalar masses reduce to : 



{mlp) = [m 2 3/2 + V - e K F™F n d m d n ln(K Q )] S, 



tt/3 



(288) 



where we have used (I267p . Using ()276l) in (|288j) . we obtain the following expression for the scalar mass 
squared 



(m 2 a ) = V + (m 2 3/2 ) [1 



A 2 P 
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where we have defined another quantity: 



(290) 
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As in the case of the trilinears, only the overall scale of the scalars can be estimated, not the individual 
masses of different flavors of squarks and sleptons. Once the cosmological constant is made small by 
imposing the constraint (|207p . the scalars are given by 



(m a ) = (m 3/2 ) 



2 

[1 " (13 5 9 2 ^3^ 6 3) ± £ V 2 n sin2 ( 27r ^) + 12 & sin ( 4 ^?) - 21 & *a(2*«?)}] 



« K /2 ) [1 - £ ^ ^ sin 2 (27r0f ) + I 2 C sin(4vr^) - 21 ft sin(27r0f )}] 

j 

~ ml /2 . (291) 

Thus, to a high degree of accuracy, in the IIA limit, the scalar masses for de Sitter vacua are flavor 
universal as well as flavor diagonal and independent of the details of the matter Kahler metric described 
by parameters 9{. Moreover, to a very good approximation, they are equal to the gravitino mass. A 
natural expectation away from the IIA limit is that the squark and slepton masses are always of order 
m 3/ / 2 - Since m 3 / 2 1S °f several TeV, the scalars are quite heavy, naturally suppressing flavor changing 
neutral currents (FCNCs). 



C. Subleading Corrections to the Condition for small Cosmological Constant and Effects on 

Phenomenology 

In this subsection we would like to give a rough estimate of how the amount of tuning of the cosmo- 
logical constant might affect the values of the soft parameters. In fact, the constraint (j207H which sets 
the cosmological constant to zero in the leading order still results in a very large value of the cosmological 
constant Vo ~ 0.01 x m^m 2 , once the subleading terms are taken into account. One can also do exact 
numerical computations for manifolds with small number of moduli (say two). 

Taking into account the terms in the subleading order as y — > 0, the potential at the minimum with 
respect to the moduli, as a function of the meson vev </>o is given by 
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Hence, vanishing of the cosmological constant corresponds to the vanishing of the combination 
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Recall that by imposing this condition in the leading order as y — * we had demonstrated that the 
combination P In ( J is fixed by 
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Thus, the leading order condition on Pin ^^^j obtained from (|294p is given by 

'A 1 Q\\ {0) 28 (Q-P) 



Pin 



AoP 



3 (Q-P) -8 



(295) 
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In this case, the meson vev 4>q in the leading order is fixed at the value 



The subleading corrections to Pin ^^^j can be found iteratively if we plug (1296P into the subleading 
term ~ y 2 in (|293p to obtain 



. 2aa 1 / abiay\ 2 f ace \ \ , 9 { aa\ 2 , _. 

Again, by setting the discriminant of the biquadratic polynomial in (|297p to zero and using 

^h+^^mKq-pW')' (298) 



we obtain the following condition on P In ( j 



28 7 / 4 4 \ . . 

+ 7 7—7 77 h- TTTTnV- ^ + 7^- ~ =»• (299) 



^ ~ P Pln(^) (pm(^)) 2 V (^) (0) (Q-^) 2 (^) (0) (Q-^ 
To compute the first subleading order correction to Pin (^^p^j we express 

where the first term in (|300p corresponds to the leading order expression in (|295p and is the subleading 
order correction. Hence, plugging ()300p into (|299f) and keeping the terms linear in after some algebra 
we obtain 

(2+(^) (0) (Q-P))(2-3(^) {0) (Q-P)) 

<5« = ^ ^ ^ . (301) 

4(<^(Q-P) 2 

Therefore, including the first subleading order, the condition on the cosmological constant to vanish 
results in the following constraint 

(A lQ \ 28(Q-P) (2 + (^) (0) (Q-P))(2-3(^^(Q-P)) 
Phl \mP) ~ 3(Q-P)-8 + 4(^)^(Q-P) 2 • (3 ° 2) 

In particular, for the case when Q — P = 3 we obtain 

P In (^j^J = 84 - 0.9977 83.002 . (303) 

which yields a ~ 1% correction to the leading order. We have confirmed that for a case when the 
compactification manifold has two moduli, numerically imposing the constraint that the cosmological 
constant has the observed value implies 

P In (Mj « 82.9958 , (304) 



68 



which agrees with (|303p to a very high degree of accuracy. 

One would now like to estimate the effects of tuning the cosmological constant to its observed value 

on phenomenological quantities. The quantity most sensitive to such corrections is the gravitino mass, 

p 

since it is proportional to ^^^^ Q P and can therefore change by a factor of order one. Of course, this 
hardly affects the distributions of scales of m 3/ / 2 and the emergence of the TeV scale peak remains very 
robust. Moving on to gaugino masses, recall that the tree-level gaugino mass is given by 

which for Q — P = 3 and the leading order constraint ^Pln ^^^p!) = ^ resulted (up to an overall 
phase) in 

(M 1/2 ) (0) » -0.0240 x m 3/2 . (306) 

Including the first subleading correction in (|303p . the tree level gaugino mass (up to an overall phase) is 
given by 



M 1/2 « -- 1 + — + — x m 3/2 = -0.0243 x m 3/2 , (307) 




resulting in a ~ 1% correction to the leading order expression in (|306|) . For a case when the compactifi- 
cation manifold has two moduli, the numerically obtained result is given by 

M 1/2 ~ -0.0242 x m 3/2 , (308) 

again confirming the high accuracy of the approximate result in (130T[) . Recall that when the leading 
order constraint Pin ^4^) ~ ^ ^ s satisfied the value of the cosmological constant is 

V ~ 0.01 m\ l2 ml ~ [O (10 10 - 10 11 ) GeV) 4 . (309) 

for m 3 / 2 ~ O (10 — 100) TeV. Thus, while the cosmological constant changes by many orders of mag- 
nitute as it is tuned to its observed value, the tuning has a very small effect on the tree-level gaugino 
mass. However, due to the cancellation between the tree-level and the anomaly mediated contributions, 
the correction to the tree-level gaugino mass computed above may be important and therefore, has been 
taken into account. ^ 

Finally, the ~ 1% correction to the leading order ^Pln f^p^ J due to the tuning of the cosmolog- 
ical constant has almost no effect on the anomaly mediated gaugino masses, trilinears and the scalars. 
This is because the terms proportional to l/Pln^^^ and 1/ ^Pln^^p^ are subleading. These 
considerations indicate that a leading order tuning of the cosmological constant is enough to allow reliable 
particle physics phenomenology. 



D. Radiative Electroweak Symmetry Breaking (REWSB) 



It is very important to check whether the soft supersymmetry breaking parameters in these vacua 
naturally give rise to radiative electroweak symmetry breaking (REWSB) at low scales. In order to check 
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that, one has to first RG evolve the scalar higgs mass parameters rrijj^ and rn/jj from the high scale to 
low scales. Then one has to check whether for a given tan /3, there exists a value of fi which satisfies the 
EWSB conditions. At the one-loop level, we find that EWSB occurs quite generically in the parameter 
space. This can be understood as follows. The gaugino mass contributions to the RGE equation for m 2 H 
push the value of up while the top Yukawa coupling, third generation squark masses and the top 
trilinear pull it down. The suppression of the gaugino mass relative to the gravitino mass causes it to 
have a negligible effect on the RGE evolution of m/jj . On the other hand, the masses of squarks and 
A-terms are both of 0(777,3/2), which guarantees that rn 2 H ^ is negative at the low scale. Typically, m 2 H ^ is 
proportional to — fTtg; 2 , up to a factor less than one depending on tan/3. Thus, the EWSB condition can 
be easily satisfied with a /j parameter also of the order 777,3/2- Note that large ^4-terms (of 0(7773/2)) are 
crucial for obtaining EWSB. Having large squark masses and small ^4-terms cannot guarantee EWSB, 
as is known from the focus point region in mSUGRA. Also, one has to ensure that the third generation 
squarks have positive squared masses, which we have checked. From the point of view of low scale effective 
theory, EWSB appears to be fine-tuned, however a better understanding of the underlying microscopic 
theory may help justify the choice of parameters. We will report a detailed analysis of these issues in 
Q 

E. The [i and B/j, problem 

We will not have much to say about the /j terms here, leaving a detailed phenomenological study for 
our future work [13] . We will however, take this opportunity to highlight the main theoretical issues. 
The normalized \x and Bfj, parameters are : 

H = (^e^V + m 3 / 2 Z- e ^ 2 F^Z)(^^J- 1/2 
W* 

Bfi = {K Hu K Hd y l/2 {^-e K l 2 ^(e K l 2 F m [K m + 8 m In//] - m 3/2 ) + (2m* /2 + V Q )Z- m^F^d^Z 

+m 3/2 F m [d m Z - Z8 m \og{k Hu k Hd )) - F™F n [drnd n Z - drnZd n \og(K Hu K Hi )]} (310) 

We see from above that the value of the physical /i and Bfj, parameters depend crucially on many 
of the microscopic details eg. if the theory gives rise to a non-zero superpotential // parameter and/or 
if a non-zero bilinear coefficient Z is present in the Kahler potential for the Higgs fields. From section 
IVIIIDl we see that one requires a /j, term of 0(m 3 / 2 ) to get consistent radiative EWSB. This is possible 
for eg. when one has a vanishing fjf parameter and an 0(1) higgs bilinear coefficient Z, among other 
possibilities. 



F. Dark Matter 

For dS vacua with a small cosmological constant, M 2 << Mi at low scale. In addition, since \i should 
be of 0(7773/2) f° r consistent EWSB as seen in section fVIII D( both M2 and M\ are much less than \i. 
Hence, the LSP is wino-like. As was discussed in section IVIII CI the tuning of the cosmological constant 
has little effect on the gaugino masses, thereby preserving the gaugino mass hierarchy. It is well known 
that winos coannihilate quite efficiently as the universe cools down. Since the wino masses in these vacua 
are 0(100) GeV, the corresponding thermal relic density after they freeze out is very small. However 
there could be non-thermal contributions to the dark matter as well, e.g. the decay of moduli fields into 
the LSP after the LSP freezes out. In addition, one should remember that the above result for a wino LSP 
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is obtained after imposing the requirement of a small cosmological constant. It would be interesting to 
analyze the more general case where the results may change. We leave a full analysis of these possibilities 
for the future 11311 . 



G. Correlations 

As we have seen, the parameters of the MSSM depend on the "microscopic constants" determined by 
a given G2 manifold and can be explicitly calculated in principle. Therefore, the parameters obtained 
are correlated with each other in general. For instance we saw that the gaugino and gravitino masses are 
related. By scanning over the allowed values of the "microscopic constants", by scanning the space of 
G<i manifolds, one obtains a particular subspace of the parameter space of the MSSM at the unification 
scale. For a given spectrum and gauge group, the RG evolution of these parameters to low scales can also 
be determined unambiguously, leading to correlations in soft parameters at the low scale Finally, these 
correlations in the soft parameters will lead to correlations in the space of actual observables (for eg, the 
LHC signature space) as well. In other words, the predictions of these vacua will only occupy a finite 
region of the observable signature space at say the LHC. Since two different theoretical constructions will 
have different correlations in general, this will in turn lead to different patterns of signatures at the LHC, 
allowing us to distinguish among different classes of string/M theory vacua (at least in principle). These 



issues, in particular the systematics of the distinguishing procedure have been explained in detail in [331 ] . 



H. Signatures at the LHC 

The subject of predicting signatures at the LHC for a given class of string vacua requires considerable 
analysis. Here, we will make some preliminary comments, with a detailed analysis to appear in [l3| . 

The scale of soft parameters is determined by the gravitino mass. We saw in section IVIIDI that 
requiring the cosmological constant to be very small by imposition of a constraint equation (eqn. (I19ip ) 
fixes the overall scale of superpartner masses to be of 0(1 — 100)TeV. Once the overall scale is fixed, the 
pattern of soft parameters at M un if is crucial in determining the signatures at the LHC. As explained 
earlier, these M theory vacua give rise to a specific pattern of soft parameters at M un if. We find non- 
universal gaugino masses which are suppressed relative to the gravitino mass. We furthermore expect 
that the scalar masses and trilinears are of the same order as the gravitino. The \i and B\x parameters 
are not yet understood. For phenomenological analysis however, we may fix them by imposing consistent 
EWSB. 

One can get a sense of the broad pattern of signatures at the LHC from the pattern of soft parameters. 
Since gaugino masses are suppressed and the fact that the anomaly contribution to the gluino mass 
parameter approximately cancels the tree-level contribution, one would generically get comparatively 
light gluinos in these constructions, much lighter than the scalars, which would give rise to a large 
number of events for many signatures, in particular many events with same-sign dileptons and trileptons 
in excess of the SM and many events with large missing energy, even for a modest luminosity of 10 fb^ 1 . 
Since the gauginos are lighter than the squarks and sleptons, gluino pair production is likely to be the 
dominant production mechanism. The LSP will be a neutralino for the same reason assuming R parity 
is conserved. 

It is also possible to distinguish the class of vacua obtained above from those obtained in Type IIB 
compactifications, by the pattern of signatures at the LHC. For the large volume type IIB vacua, the 



scalars are lighter than the gluino [34;] while for the KKLT type IIB vacua, the scalars are comparable to 



the gluino 35[. This implies that squark-gluino production and squark pair production are the dominant 



production mechanisms at the LHC. Since the LHC is a pp collider, up- type squarks are preferentially 
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produced from t-channel valence u-quark annihilation if they are reasonably light, leading to a charge 
asymmetry which is preserved in cascade decays all the way to the final state with leptons. On the 
other hand, for the class of M theory vacua described here, since gluino pair production is the dominant 
mechanism and the decays of the gluino are charge symmetric (it is a Majorana particle), the M theory 
vacua predict a much smaller charge asymmetry in the number of events with one or two leptons and > 
2 jets compared to the Type IIB vacua. 



I. The Moduli and Gravitino Problems 

The cosmological moduli and gravitino problems can exist if moduli and gravitino masses are too 
light in gravity mediated SUSY breaking theories. Naively, after the end of inflation, the moduli fields 
coherently oscillate dominating the energy density of the universe. Since the interactions of the moduli 
are suppressed by the Planck scale (m p ), their decay rates are small perhaps leading to the onset of a 
radiation dominated universe at very low temperature (T R ~ 0(1(T 3 ) MeV for moduli of 0(100 GeV- 
5TeV)), compared to what is required for successful BBN. 

To check if the moduli and gravitino problem can be resolved in these M theory compactifications, 
one has to first compute the masses of the moduli. After doing this, we will use the results to discuss the 
moduli and gravitino problems. 

The geometric moduli Sj appear in the lagrangian with a kinetic term given by 

N 

^j^d^Si, (311) 
i=i Si 

which is non canonical. The canonically normalized moduli Xi are 



\ ^In.s,. (312) 



The complete mass matrix including the mixed meson-moduli entries is given by: 

, 2 \ 2 zv 2 (aj Oj)^ d 2 V , 2 \ /2aj v d 2 V , 2 \ 9 2 V 

\ m x)ij = q AT. AT. Ae.Ae. ' l m )id>„ = V ~5~W~a7^Z ' V™ ' </»o 4>a = ' ^ 33 ' ) 



where we took into account the fact that = 0, |^ = at the extremum and that Vi = v for all 
% = 1, N. A fairly straightforward but rather tedious computation yields the following structure of the 
mass matrix 

( m2 x)ij = (( a * a ^ K i + s a K *) x m 3/ 2 » ( 314 ) 
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where, in the large v approximation, K\, K 2 , K3, K4 are given by 
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Using condition P In y-^pj = 84 to tune the cosmological constant together with Q — P = 3, we obtain 
from (HDU) 



K , ^ SGO 1 G ( |] , K 2 ~ 3.83 , K 3 « 6815 , K 4 « 540 f ^ 



(316) 



We can diagonalize the matrix ()314p in two steps. We first construct a set of orthogonal (unnormalized) 
N + 1 vectors given by 
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Next, we normalize the above vectors and construct an orthogonal {N + 1) x (N + 1) matrix 



x\ x 2 
\x\\ ' \x 2 \ ' 



XN+l 

\xn+x\ 



(317) 



By applying the orthogonal transformation 1Z, the mass matrix in (|314p is converted into a block-diagonal 
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form given by: 



(K 2 
K 2 










V 








K 2 
\K X + K 2 







1*3 







1*3 
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Hence, the first N — 1 eigenvalues are given by 

Xi = K 2 for i = 1,...,N -1. 



(318) 



Diagonalizing the 2x2 block is then trivial and boils down to solving a simple quadratic equation. Since 
for typical values of Q and P we have 



the corresponding eigenvalues are given by 



K 3 ^ , K 2 

— - < 0.07 , and -f- < 1 , 



#1 



(319) 



(320) 



Because of the significant cancellation due to the minus sign in Xn+i we have to beyond the leading order 
approximation in (|316p when computing K\, K 3 and K±. For dS vacua with a nearly zero cosmlogical 
constant we have verified numerically that the corresponding eigenvalue is always positive, confirming 
that we have obtained a stable minimum. The N — 1 degenerate light states have masses given by 



K 2 x m 3/2 « 1.96 x m 3/ /2 , k 



1,...,N 



1 



(321) 



independent of any parameters of the model. For dS vacua with a very small cosmological constant, 
numerical computations for the choice P = 27, Q = 30 give the following masses for the remaining two 
states 



Mat « 600 x m 3 / 2 , Mn + i 2.82 x m 3 / 2 



(322) 



Note that all the masses are independent of the number of moduli N as well as the rationals a%. 

Since the gravitino mass distribution peaks at 0(100) TeV, which is also in the phenomenologically 
relevant range, and the light moduli are roughly twice as heavy compared to the gravitino, the moduli 
masses are heavy enough to be consistent with BBN constraints. A detailed study of this issue will 
appear in 



IX. SUMMARY AND CONCLUSIONS 



A major goal of string/M theory is to find solutions that incorporate the Standard Model, important 
clues to physics beyond the SM such as gauge coupling unification, and in addition explain phenomena 
the SM cannot explain. The most important unsolved problem is explaining the value of the weak scale 
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(at or below about a TeV) purely in terms of the Planck Scale of O(10 19 GeV) or some other fundamental 
scale, the Hierarchy Problem. One obstacle to solving the hierarchy problem is that all the moduli that 
characterize the string theory vacuum must be stabilized for a meaningful solution. 

In fluxless M theory vacua the entire effective potential is generated by non-perturbative effects and 
depends upon all the moduli. In this paper we have studied this potential in detail for a particular form 
of the Kahler potential when the non-perturbative effects are dominated by strong gauge dynamics in 
the hidden sector and when such vacua are amenable to the supergravity approximation. In the simplest 
case, we studied GVmanifolds giving rise to two hidden sectors. The resulting scalar potential has AdS 
vacua - most of them with broken supersymmetry and one supersymmetric one. Then we studied the 
cases in which there was also charged matter in the hidden sector under the plausible assumption that 
the matter Kahler potential has weak moduli dependence. In these cases the potential receives positive 
contributions from non-vanishing F-term vevs for the hidden sector matter leading to a unique de Sitter 
minimum. The de Sitter minimum is obtained without adding any "uplifting" terms, such as those 
coming from antibranes, which explicitly break supersymmetry. With the form of the Kahler potential 
given by ([I]), we have explicitly shown that all moduli are stabilized by the potential generated by strong 
dynamics in the hidden sector. 

In the de Sitter minimum we computed 7713/2 and found that a significant fraction of solutions have 
777,3/2 in the TeV region, even though the Planck scale is the only dimensionful parameter in the theory. 
The suppression of 7773/2 is due to the old idea of dimensional transmutation which has also been used 
in heterotic string theory. The eleven dimensional M theory scale turns out to be slightly above the 
gauge unification scale but below the Planck scale. The absence of fluxes is significant for simultaneously 
having 777,3/2 ~ TeV naturally and M\\ not far below the Planck scale. 

The problem of why the cosmological constant is not large is of course not solved by this approach. 
We do however understand to a certain extent what properties of £r2-manifolds are required in order 
to solve it and we suggest that one can set the value of the potential at the minimum to zero at tree 
level and proceed to do phenomenology with the superpartners whose masses are described by the softly 
broken Lagrangian. One particularly nice feature is that we are able to explicitly demonstrate that the 
soft-breaking terms are not sensitive to the precise value of the potential at the minimum. A tree-level 
tuning of the cosmological constant is enough to determine the phenomenology. 

When we set the value of the potential at the minimum to zero at tree level a surprising result occurs. 
Doing so gives a non-trivial condition on the solutions. When this condition is imposed on 7773/2, for 
generic G2 manifolds it turns out that the resulting values of 7773/2 are all in the TeV region. Thus we do 
not have to independently set Vq to zero and set 7773/2 to the TeV region as has been required in previous 
approaches. 

A more detailed study of the phenomenology of these vacua, particularly for LHC and for dark matter, 
is underway and will be reported in a future paper. In the present paper we presented the relevant soft- 
breaking Lagrangian parameters and mentioned a few broad and generic features of the phenomenology, 
for both our generic solutions and for the case where Vo is set to zero at tree level. We presented a 
standard supergravity calculation of the soft breaking Lagrangian parameters, and found that the scalar 
masses m a , and also the trilinears, are approximately equal to 7773/2, to the extent that our assumptions 
about the matter Kahler potential are valid. Remarkably, the tree-level gaugino masses are suppressed 
by a factor of 0(10 — 100). This suppression is present for all GVmanifolds giving the de Sitter minimum. 
For calculating the tree level gaugino masses the matter Kahler potential does not enter, so the obtained 
values at tree level are reliable. Because the gaugino suppression is large, the anomaly mediated mass 
contributions are comparable to the tree level ones, and significant cancellations can occur. Gluinos are 
generically quite light, and should be produced copiously at LHC and perhaps even at the Tevatron - this 
is an unavoidable prediction of our approach. We have also checked that radiative electroweak symmetry 
breaking occurs over a large part of the space of G2 manifolds, and that the lightest neutralino is a good 
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dark matter candidate. It will be exciting to pursue a number of additional phenomenological issues in 
our approach, including inflation, baryogenesis, flavor and CP-violation physics, Yukawa couplings and 
neutrino masses. 

The approach we describe here apparently offers a framework that can simultaneously address many 
important questions, from formal ones to cosmological ones to phenomenological ones (apart from the 
cosmological constant problem, which might be solved in a different way). Clearly, however much work 
remains to be done. In particular, a much deeper understanding of ^-manifolds is required to understand 
better some of the assumptions we made about the Kahler potential of these vacua. 
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APPENDIX A: KAHLER METRIC FOR VISIBLE CHIRAL MATTER IN M THEORY 

As stated in section IVIIIBl we will generalize the result obtained for the Kahler metric for visible 
sector chiral matter fields in toroidal orientifold constructions in IIA [2^ | to that in M theory. The 
result obtained for the Kahler metric for the (twisted) chiral matter fields (cj> a ) in the supergravity limit 
(a' -> 0) in [28] is: 



TT ( HL~ ^ x 1/2 



K % = 111 ~ r(0«V ) S ^ 



i=i 

tanfrfl?) = - ; lP = Vi + iT4 (Al) 

U{ + qf/pf 

where i = 1, 2, 3 denote the number of moduli in the specific IIA example, {p, q} are integers and U l are 
the complex structure moduli in type IIA in the geometrical basis. As mentioned before, we will restrict 
to factorized rectangular tori with commuting magnetic fluxes, for which U{ = 0. Then 

tan(7<) = ^ (A2) 

The first step towards the generalization is to identify the G2 moduli in terms of IIA toroidal moduli 
by imposing consistency between results of IIA and M theory. The consistency check is the formula for 
k 2 - the physically measured four dimensional gravitational coupling. We have: 

- 2 = voW MTheory W- 

K 2 = ^ut\ > HA String Theory (eqn. 18.2.2 of [36]). (A3) 
Vol(A6 ) 



76 



where X7 and Xq are the volumes of the internal 7-manifold and 6-manifold (in IIA) respectively. The 
M theory gravitational coupling and the string theory gravitational coupling k 2 can be written in 
terms of the string coupling in IIA {g s = e^o) and a . 



k 2 



K 2 



1 ,~ ,7 /4 



-(2tt) V (eqn. 13.3.24 of [36]) 



(A4) 



Also, the volumes of Xj and Xq (for a IIA toroidal orientifold) can be written as: 

N 

(sA a '- 

47T 



;9 

Vol(X 7 ) = V x l 7 M where V x = His^; -M- = k 2 u [18] 



i=l 



Vol(X 6 ) = (2vr4 j )(2vr4 j )(2vr4 , )(2vr J Rr ) )(27r J R^)(2^4 



Using (|A4p . we get 



Im 



which gives rise to the following expression for k in M theory: 



/ 2/3 
vra g s ' 



V x 



(A5) 
(A6) 

(A7) 



In IIA String theory, the definition of the IIA moduli Tj and U{ in terms of the geometry of the torus 
is given below. We will stick to the case of a factorized T 6 , rectangular tori, and commuting magnetic 
fluxes (in the type IIB dual) for simplicity, in which case only the imaginary parts of the moduli are 
important: 



Im(TY 



R 



R\ 



(0 



im{uy = 



R\ K 2 

a' 



1,2,3. 



(A8) 



(i) (i) th 

where R\ , R% are the radii of the i IIA torus along the x and y axes respectively. 



Now, from (|A3]> . (IA4I) . (jA5J) and (fASjl we get: 



Combining (|A7jl and (IA9I) gives: 



na'g 2 



(A9) 



9 r 



(A10) 



The above formula (eqn. (|A10|) ) is quite general and should always hold 9 , since it has been derived by 
requiring consistency between formulas for the physically measured gravitational coupling constant. To 
identify the individual G2 moduli however, is harder and model-dependent. In the next subsection, we 
will do the mapping for the case of a simple toroidal G2 orbifold (T 7 /Z$) considered in [37J, where it can 
be shown that eqn (|A1Q[) is satisfied. 



within the limits of the IIA setup considered. 
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1. Particular Case 



In 



371 ]. the definitions of the G2 moduli (eqn 2.4) and the Kahler potential (eqn 2.10) are not given 



in a dimensionless form. Therefore, we will make them dimensionless as is done in [181 ]. So, we have: 



- r = ^—- if = -31n( ^ )+ constant (All) 



The volume of the manifold and the moduli in this compactification are explicitly given as 371 ]: 

7 



Vol(^7) — |^ a} — R-LR2R7] o 2 — RiR'iRa] — R1R4R5] 
1=1 

= R2R3R5', a = R2R4RQ] 0^ = R3R4R7] = R^RqRj • (A12) 



From (|A12|) . we can write Vol(X 7 ) = {a 1 a 2 o? a^a 5 a Q a 7 ) 1 ^ , which implies that Vx = nj=i( s «) 1 ^ 3 - There- 
fore, in the notation of [lgj, = 1/3; i = 1, 2, 7. We will identify the M theory circle radius as R7, 
the remaining six radii can just be identified as the x and y radii of the three tori in IIA: 

Rx = (27r)i?5 1) ; R 2 = {2tt)R { 2 1] ; R 3 = (2vr) J Rf ) ; R 4 = (2tt)^ 2) 

R 5 = (2vr)i?f ) ; i? 6 = (2tt)R^; R 7 = (2ir)g s a' 1/2 . (A13) 



With these identifications, and using (|A6|) . we can write the individual G2 moduli in terms of the IIA 
moduli: 

1 /ri%«c/f t/f V /2 . _ 1 N 2) uPuPu? V /2 . 



52 5, ^ T 2 (1) T 2 (2) J 5 53 9s \ T 2 (1) T 2 (3) ; ' 

_ 1 / jg^rjf^g) y/ 2 ! 1/2 

si V ???? / ^ ' ' 

Therefore, we get: 

7 (l) ( 2 ), r (3) 

Vx = II(^) 1/3 = 2 4% 2 ( Al4 ) 

i=l ^ 

which is the same as (lAlOj) . So, for the particular case, this suggests the following generalization: 

T(l-0f^ 1/2 



K a = n 

i=l,6,7 



r(e?) 



tan(7T0?) = z = 1,6,7 (A15) 



2. General Case 

For more general G2 manifolds with many moduli, the precise map of the individual moduli is not 
completely clear. However, it seems plausible that the complex structure moduli appearing in the Kahler 
metric in IIA map to a subset of the G2 moduli in a similar way, as in the particular case. 
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Therefore, we use the following expression for the Kahler metric for visible chiral matter fields in M 
theory: 



n 



r(i 



Yip?) 



tan(-7r6 



i=l 



1/2 



l,2,..,p<N(=N) 



cf (sj) ; cf = constant; I = rational number of 0(1). 



(A16) 



The derivatives of the Kahler metric w.r.t the moduli are very important as they appear in the 
soft scalar masses, the trilinears as well as the anomaly mediated gaugino masses, as seen from section 
IVIIIB1 The first derivatives appear in the trilinears and anomaly mediated gaugino masses while the 
second derivatives appear in the scalar masses. For the metric in (IA16I) . these can be written as follows: 



d n ln(K Q )=^ 



drnd n ln{K a 



era 



901 

dz n 



1 d 

2d9°, 



ln( 



rfi 



dz n dz 



4(C) = (A17) 



The functions tfj® an< ^ ^nn depend on the angular variable 9", which in turn depend on the moduli. The 
first and second derivatives of 9® with respect to z n are given by: 



1 



89% 



I 



4iri s n 



del 

dz n 

dz n dz n ' 4 d(s n ) 2 16vr (s n ) 2 



2i ds n 

i d 2 e° 



sm(27r0£) 



(A18) 



/ 



[Zsin(47r^) -2sin(27r0°; 



The dependence of the soft parameters in section IVIIIBI on 9% is extremely simple. Instead of re- 
expressing the dependence on 9" in terms of the moduli, it is much more convenient to retain the 
dependence on as 9" E [0, 1) and the variation of relevant functions with 9® in the allowed 
range can be plotted easily. In particular, the function F(9%) = ^{ip" srn (27r#n)} appears in the ex- 
pression for the anomaly mediated gaugino masses (eqn. (|279|) ) and the trilinears (eqn. (|287|) ). the 

function G(#") = In ^ ^ - ^ appears in the expression for the trilinears (eqn. (|287p ) and the function 

H{ d n) = ^{l 2 ^m sin 2 (2vr6»°) + I 2 ip% sin(4vr^) - 2/^n sin(2vr^)} appears in the expression for the 
scalars (eqn. (|291fl ). The variation of these functions with 9" is quite mild as seen from Figure lA2l 




FIG. 21: Left: F as a function of 9%. Middle: G as a function of 6%. Right: H as a function of 9* (for 1=1). From J2g], 

9 n e [0, l). 



Since the functions F, G and H vary very mildly with 9" and are 0(1) in the whole range, it is 
reasonable to replace the above functions by 0(1) numbers, as is done in section IVIIIBI This is justified 
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since we are only interested in estimating the rough overall scales of various soft parameters. 
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